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Abstract: We study the mini-superspace quantization of spatially homogeneous (Bianchi) cosmological uni- 
verses sourced by a Dirac spinor field. The quantization of the homogeneous spinor leads to a finite-dimensional 
fermionic Hilbert space and thereby to a multi-component Wheeler-DeWitt equation whose main features are: 
(i) the presence of spin-dependent Morse-type potentials, and (ii) the appearance of a q-number squared-mass 
term, which is of order 0(h 2 ), and which is affected by ordering ambiguities. We give the exact quantum 
solution of the Bianchi type-II system (which contains both scattering states and bound states), and discuss 
the main qualitative features of the quantum dynamics of the (classically chaotic) Bianchi type-IX system. We 
compare the exact quantum dynamics of fermionic cosmological billiards to previous works that described the 
' spinor field as being either classical or Grassmann-valued. 

^ ■ 1 Introduction 

The main aim of this work is to clarify the physical structure of the coupled quantum Einstein-Dirac system 
within a minisuperspace cosmological setting, and, in particular, its dynamics in the vicinity of a cosmological 
singularity. Let us recall that the seminal work of Belinsky, Khalatnikov and Lifshitz (BKL) |TJ (see also 
Misner [2]) has brought into light the chaotically " oscillatory" behaviour, near the singularity, of the diagonal 
components of the metric both in homogeneous Bianchi IX cosmological models, and in the "general (classical) 
solution" of Einstein equations. Recently, this chaotic BKL behaviour acquired a new significance through the 
discovery of its unexpected link with infinite-dimensional Kac- Moody algebras [21 HI [5] (for reviews, see [Hl[7])- 
The presence of such a chaotic behaviour depends both on the spacetime dimension, D, and on the matter 
content of the considered cosmological model. For instance, pure gravity (i.e. the vacuum Einstein equations) 
has a chaotic, BKL-like behaviour in spacetime dimension D < 10, but a monotonic, Kasner-like behaviour in 
spacetime dimensions D > 11 [§]. 

The effect, near a singularity, of a general bosonic matter content (scalar field and p— forms), in any spacetime 
dimension D, has been studied in Ref. [6]. The main result is that almost all the bosonic degrees of freedom 
"freeze" near the singularity (i.e. admit limits, at each spatial point, as t — > 0), except for the diagonal part 
of the spatial metric, parametrized as g aa = exp(— 2j3 a ) (together with any scalar field tp = /3q, if present). 
The dynamics of the /3 a 's can be represented as that of a "ball" (of position f3 a ) moving on a Lorentzian (or 
hyperbolic) billiard. The latter billiard motion can then either be chaotic, or non-chaotic (i.e. ultimately 
monotonic), depending on the geometry of the " billiard table", which is a polyhedron domain in hyperbolic 
space. The quantum mini-superspace versions of some of these bosonic cosmological billiards have been studied 
in several works [51 UD1 ITU IT2] 

While the effect of bosonic matter on the dynamics near a singularity is well understood, the effect of 
fermionic matter is less well understood. This difference in understanding has both technical and physical roots. 
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From the technical point of view, the Hamiltonian description of spinor fields, coupled to gravity, is rather subtle 
and complex. The Hamiltonian description of a Dirac (spin i) field, coupled to gravity (notably in Bianchi 
spacetimes), has been clarified in several papers, see Refs. [HI HH [15] [HI [TT1 [HI HU [20]. Here, we shall not 
consider the case of the gravity coupling to a Rarita-Schwinger (spin |) field, i.e. the case of "supersymmetric 
(quantum) cosmology" , though we view the Einstcin-Dirac system as a toy model for the Einstein-Rarita- 
Schwinger case. See, e.g., [21] for a entry in the literature on supersymmetric quantum cosmology. On the 
other hand, the physical meaning of having spinorial sources in cosmology has remained somewhat obscure. 
Various (conflicting) approaches to this issue have been assumed in the literature. We shall distinguish three 
different approaches to the treatment of spinorial fields, say ip: (i) bilinearfQ in ip (in the source terms, , for 
gravity) are replaced by real numbers (or "c-numbers" ) ; (ii) the ip variables are treated as Grassmann- valued 
(or "G-numbers" ) ; or (iii) the ip's are treated as quantum operators (or "q-numbers"). 

When treating the ip's as G-numbers, they do not affect the BKL chaos (if the latter is present in absence of 
fermions). Indeed, T M " oc ip 2 only modifies the "soul" 32+54 + - ■ • of the metric 3, i.e. the part of the Grassmann 
expansion 3 = go + 32 + 34 ■ ■ • of g which contains an even (non zero) number of Grassmann generators. By 
contrast, the "body" go of g (i.e. the part of g which does not contain any Grassmann generator) is unaffected 
by the ip's, and is entirely determined by the bosonic sector of the theory. This approach has been recently 
used 22 to discuss how the chaotic behaviour of the body go of the metric induces a corresponding chaotic 
dynamics in the (Grassmann-level-one) fermions ip — ip\ + tp3 + ■ ■ ■ • 

By contrast, if one replaces the ~ ip 2 source terms by c-numbers (that react back on the body 30 of the 
metric), they can drastically modify the asymptotic behaviour near the singularity. Indeed, it has been argued 
by Belinsky and Khalatnikov [53] that the presence of a Dirac field would, if so treated, ultimately destroy 
the BKL chaos. This result has been confirmed, and streamlined, by the Hamiltonian treatment of [24] which 
showed that the ~ ip 2 source terms modify the billiard dynamics of the (logarithmic) cosmological scale factors 
[3 a by adding a (positive) "squared-mass term" /z 2 to the effective (Lorentzian-signature) Hamiltonian describing 
the dynamics of f3 a . Indeed, such a squared-mass term (in Lorentzian j3 space) slows down the motion of the (3 
particle and ultimately prevent collisions on the cushions of the billiard table. 

In the present work we shall treat the fermions as q-numbers and face the problem of discussing the meaning 
of the back reaction of fermions within a quantum framework. This is a notoriously difficult problem, but 
we shall be able to get answers by restricting ourselves to the mini-superspace framework of homogeneous 
cosmological spacetimes. In that framework, the quantization of the spinor ip is equivalent to considering that 
the wavefunction $ of the universe is a multi-component object, say $ CT , where the discrete index a labels the 
finite number of independent states allowed by the anticommutation relations of ip. In particular, as we shall see 
in detail below, when ip is a Majorana spinor the discrete index a takes 4 values and can be identified with the 
spinor index of an auxiliary 50(4) Clifford algebra defined by the quantum operators T a = 2g*ip a , a — 1,2, 3, 4. 
The appearance in the wavefunction of such a discrete spinor-type label when considering zero-mode (spatially 
independent) Fermionic operators is well known from the Ramond string case |25j , and has also been used in 
some studies of quantum supersymmetric minisuperspace cosmologies, see, e.g., |26j . 

Most of our discussion will allow for a general (class A) Bianchi model, but we shall discuss the quantum 
dynamics in details only for two cases: the Bianchi IX (or 'mixmaster') model, and the (non-chaotic) Bianchi 
II model. We shall also compare the quantum solution to the various ways of discussing its classical analogs. 

Our paper is organized as follows. Section 2 discusses the classical Lagrangian formulation (with a specific 
fixing of the vielbein in terms of the metric) of an homogeneous spinor field coupled to a Bianchi metric. 
This is followed in Sec. 3 by the corresponding classical Hamiltonian formulation. The quantization of this 
system, and the discussion of the quantum dynamics of the Bianchi- IX and Bianchi-II Einstein-Dirac systems 
is presented in Sec. 4. The comparison of the q-number and G-number approaches is discussed in Sec. 5, 
while our main conclusions are presented in Sec. 6. Finally, three appendices present complementary material, 
namely: Appendix A : alternative Hamiltonian approach; Appendix B: alternative fixing of the dreibein; and 
Appendix C: classical dynamics of the Bianchi-II case. 

iSome authors even replace the -0 variables themselves by real numbers. 
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2 Classical Lagrangian formulation of a homogeneous spinor field 
coupled to a Bianchi cosmological metric 

2.1 Choice of approach 

The coupling of a spinor field to gravity poses special problems in view of the need to use a vielbein, say h a ^ 
in addition to the metric g^ u , to describe the spinor degrees of freedom. The metric and vielbein components 
are related by 

9^ = r) & 0h%hi. (2.1) 

Here, and in the following, we shall use hatted indices to denote "flat indices" referring to a local orthonormal 
frame where r/ & a = diag(— 1, +1, +1, +1). [We restrict ourselves to a four-dimensional spacetime, and use greek 
indices to denote spacetime indices, and latin indices to denote spatial ones.] 

There are basically two different approaches to the description of the coupled gravity-spinor (or Einstein- 
Dirac) dynamics. Either: (i) the gravity degrees of freedom are only described by means of the metric com- 
ponents g^, the vielbein ones h a ^ being (locally) determined in terms of <7 M „ by some suitable gauge- fixing of 
the local Lorentz symmetry 5*0(3, 1), or (ii) one does not break the local Lorentz symmetry, and describes the 
gravity degrees of freedom by means of the redundant vielbein variables h" . The approach (i) is technically 
simpler, but depends on the choice of a specific gauge-fixing of the local 5*0(3, 1) symmetry. The approach (ii) 
does not depend (until one discusses explicit solutions) on the choice of a gauge-fixing of 50(3, 1), but involves 
more constraints than the first approach, namely the constraints linked to the local 50(3, 1) gauge symmetry. 
We shall use the approach (i) in the text, and discuss the approach (ii) in an Appendix. Note that the approach 
(i) has been advocated in Ref. [17], but that its implementation in that reference differs from the one we shall 
use (and does not take advantage of the useful automorphic potentialities of the matrix 5 discussed below). 

Having decided on the approach (i), we need to choose a specific way of fixing the local 50(3,1) gauge 
symmetry, i.e. of determining a specific vielbein /i Q , when given the metric components <7 M „ (in some coordinate 
system, or, more generally, some non-holonomic frame). Here, the (assumed) symmetry properties of Bianchi 
models come to our help. Let us recall that the geometry of homogeneous spacetimes admits the special form 

ds 2 = g^v (x) dx^ dx v 

= -N 2 {t)dt 2 + g ab (t)(N a (t)dt + T a (x))(N b (t)dt + T b (x)) (2.2) 

where the one- forms r a (x) = r° , (x)dx t only depend on the spatial coordinates used in the t = const, slices. 
In the case (considered here) of Bianchi geometries, i.e. such that the spatial slices admit a simply transitive 
Lie group G preserving the geometry, the one-forms r° can be chosen to be invariant under the group G, i.e. 
£^ h T a = 0, where £ denotes a Lie derivative and = Q(x) d/dx l , with b = 1, 2, 3, a basis of three infinitesimal 
generators of G. With a suitable choice of r a , the structure constants C c a b of (the Lie algebra of) G, which 
enter the Lie brackets of the Killing vectors, [£ a , £&] = +C C a b £c, also enter the Cartan differential of the forms 
r a , namely 

dr a = + l - C a bc T b Ar c . (2.3) 

Note that the constants C a be enter with the same sign in £ c ] and in dr a , because there are two flips of sign 
when going from the £ a 's to the r a, s: one flip between the bracket of [£, £] and the d of the co- frame dual to the 
£'s, and a second flip between the structure constants of the left action of G (defined, say, by the £'s) an d those 
of its right action (which commutes with the £'s, and corresponds to the invariant (co)-frame linked to the r's). 

Apart from the r a 's (which depend on the spatial coordinates only), the other objects entering the homoge- 
neous metric (|2.2j) only depend on the (coordinate) time, t. In addition, by using some time-dependent change 
of coordinates, one can ensure that the components of the shift vector N a (t) vanish for all times. We shall 
generally assume that this is the case, but we have left them in Eq. (|2.2|) as a reminder that, in the Hamiltonian 
formalism, the A a 's enter the action as Lagrange multipliers of the three diffcomorphism constraints H a ~ 0. 
Similarly, the "lapse" N(t) enters the action as a Lagrange multiplier for the Hamiltonian constraint T-L ra 0. 
Finally, the dynamical variables of Bianchi geometries are the six functions of time g a b(t). 
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2.2 Fixing the local Lorentz gauge symmetry 

In view of the special structure (|2.2p of the geometry it is natural to choose a vielbein co-frame 8 a — h a „ dx^ 
of the form 

6> 6 = N(t) dt , 

9 & = Y< e ~ Mt) S& b {t){T b + N b (t)dt) 1 

b 

where the matrix S a b satisfies^ 

5>" 2 ^ S & b S' a c = g bc . 

a 

In other words, the matrix S a or rather its inverse S b a, such that 

S & c S c b = 6f, S a £ S £ b ^S^, (2.7) 

transforms the quadratic form g bc into a diagonal form 

g ab S a i S b j = [diag(e- 2 ^ , e" 2 ^ , e" 2 ^ . (2.8) 

The idea behind the representation (|2.6I) is to encode the six independent components of g ab into two 
sets of three variables: (1) three "diagonal" degrees of freedom 01,^2, 03', and (2) three "off diagonal" ones, 
parametrizing the "diagonalizing" matrix S. For such a decomposition to be uniquely defined, one needs to 
restrict the structure of the (3 by 3) matrix S by six conditions. This can be achieved in several different 
ways. For instance, one could require the lower-than-diagonal elements of the matrix S a b to vanish, and its 
diagonal elements to be equal to 1; this would correspond to the so-called Iwasawa decomposition of h a b (and, 
correspondingly, of g c d = S~t h a c h b j), which is unique, and which was found to be useful in recent work on 
the hidden presence of Kac-Moody symmetries in gravity theories [5] . We shall discuss in Appendix B the use 
of this decomposition in the study of the dynamics of type II Bianchi cosmologies. However, in the case of the 
most generic (class A) Bianchi models, namely type IX and type VIII, the Iwasawa decomposition is rather 
inconvenient. 

As was emphasized by R.T. Jantzen [HI [27l [28] it is quite advantageous to use decompositions (|2.5|) . 
(12. 6|) with a matrix S a b restricted to belonging to the automorphism group, say A, of the Lie algebra 25 of G. 
Explicitly, this means that A is the group of linear transformations which leave invariant the structure constants 
C a bc . The reason why the choice of a "diagonalizing" matrix S belonging to the automorphism group A is 
advantageous is that (as we shall see explicitly below) the potential terms in the Hamiltonian can be expressed 

in terms of the /3's and of the components, say C" -r-, of the structure constants w.r.t. the intermediate, co-frame 
T a = S a b r b . These components are given by 

C" l - = S & a/ S b ' B S c ' £ C a ' b , d . (2.9) 

For a general choice of S a b the components C" would depend on the off-diagonal variables entering S a b . 

For instance, when using an Iwasawa decomposition of g ab (with an upper diagonal matrix S), the 
explicitly depend on the off-diagonal variables ^12, ^23, ^13 entering Eq. (|B.1|) . so that the potential terms in the 
Hamiltonian also depend on these off-diagonal variables. By contrast, by definition of the automorphism group 
A (as fixing © and therefore the C"s) when S £ A the components C° j- are simply equal to the original C a bc , 
and thereby do not introduce any dependence on the off-diagonal metric variables. 

When G is a simple group (and when its Dynkin diagram has no symmetries), the automorphism group A 
of © is the adjoint group of G. In the case of Bianchi type IX (where G = SU(2)), this means that A is SO(3) 
(which is the quotient SU(2)/Z2). In that case, S a b is an orthogonal matrix, and the decomposition (|2 . 6[) is 
the so-called "Gauss decomposition", corresponding to the diagonalization of the quadratic form g ab w.r.t. a 

2 We generally use Einstein's summation convention, except when there are ambiguities, as in Eqs. 1 12. 5 I t or ( 12.61 ). 
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(2.4) 
(2.5) 

(2.6) 



given Euclidean metric 8 a b- [Such a Gauss decomposition was advocated by M. Ryan [lOj [TT].] In the case 
of Bianchi type VIII, the automorphism group is 50(1,2), which means that one should use an "hyperbolic" 
generalization of the Gauss decomposition of g a b w.r.t. a given Lorentzian metric rj a b = diag(— 1, — 1, +1). The 
special metrics, 8 a b or rj a b, that enter here are simply (modulo a suitable normalization) the Cartan-Killing 
metrics fc b associated to the Lie algebra 0, say 

k ab = -\ C * ad C\ c . (2.10) 

For instance, in the usual basis for type IX we have C a be = Sabc, so that k a b = + 5 a b- 

Summarizing: In the Bianchi type IX case, we parametrize the six metric degrees of freedom contained in 
9ab{t) by means of the three diagonal variables (3i(t), /^(i), Pzit) and the three Euler angles 0i(t), 02(t), Qz{t) 
parametrizing the orthogonal metric S a b (81,62, 63) entering the Gauss decomposition (|2.6|) of the quadratic form 
9ab ("transformation to principal axes"). On the other hand, in the Bianchi type VIII case, the parametrization 
of g a b(t) by three diagonal variables (3i(t), ^(t), ^(t) and three "diagonalizing angles" 0%(t), 62(t), 03(t) should 
be done by an hyperbolic (S € SO(l, 2)) generalization of the Gauss decomposition, i.e. a transformation of g a b 
to principal axes w.r.t. a Lorentzian metric 77,^ = diag(— 1, — 1, +1). [Such a transformation is always possible. 
It can be built from the eigenvectors of the non positive-definite quadratic form r]^ w.r.t. the positive-definite 
one g a b ■} 

Then, in terms of such a parametrization, g a b *H> (j3i , fa , @3 , S a b (^1,^2,^3)) of the metric, we gauge-fix 
the local Lorentz symmetry by defining the specific vielbein 6 a (f3 a , 9 a ) = h^(f3 a , 9 a ) dx^ by means of Eqs. (|2.4p . 
(1231) . 



2.3 Lagrangian formulation of the Einstein-Dirac system 

Having uniquely determined (for types IX and VIII) a vielbein a ^ dx^ in terms of the usual metric degrees of 
freedom g a b(t), N a (t), N(t), we can now consider the general Einstein(-Hilbert)-Dirac Lagrangian density, 

L = L EH + L Dl (2.11) 

wher^l 



Leh = V^9 A R, (2-12) 

Ld = \/~ T g (fT^V^-mf^) , (2.13) 

as a functional of the metric g^, and of the spinor field * (w.r.t. the gauge-fixed vielbein h"(g„\)). We use 
gamma matrices adapted to our mostly plus signature — h H — h, namely 

7^+7^7" =277^1, (2.14) 

with an anti-hermitian 7 (satisfying (7 ) 2 = —2), and three hermitian 7 1 ,7 2 ,7 3 . We then choose 

/3:=i 7Q = -i7 5 , (2.15) 

with the involutive /3 (0 1 = +11) being used to define the usual Dirac conjugate (as defined in the mostly minus 
signature) 

*:=*t£. (2.16) 

The frame covariant derivative of the spinor entering the Dirac action Eq. (j2 . 13[) is V„ ^ = & V M where 
a is the vielbein frame dual to the vielbein co- frame h a ^ (i.e. h a ^ ^ — S a & ) , and where the world- index 
covariant derivative V^, "J is given by 

V^ = d^ + \u & ^j & ^, (2.17) 



'We use units such that 167r G = 1 = c. 
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where (denoting h afl = r) & ah^ M ) 



are the connection components of the vielbein (with last index taken as world index), and wher^l 

:= 1 (7*7* -7^7 a ) = 7 [< V ] ■ (2.18) 

We shall only consider here class A Bianchi models, i.e. models satisfying C a ac — 0, which is equivalent to 
saying that the dualization of the structure constants w.r.t. the antisymmetric lower indices, n ad :— | e bcd C a be, 
yields a symmetric tensor density. In other words 

C a bc = e bcd n ad , (2.19) 

where e a bc = £[ a bc] (with £123 = +1), and n ab = n ba . For type IX, one has n ab — S ab in the usual basis, 
while, for type VIII n ab = diag(+l, +1, — 1). [Note that the Cartan-Killing metric k a b, Eq. (|2 . 10[) . associated 
to the C's is quadratic in n ab . In type IX k a b — S a b — diag(+l, +1, +1) is numerically equal to n ab , while in 

type VIII k a b — diag(— 1, —1, +1) is of signature h, independently of whether one chooses a basis where 

n ab = diag(+l, +1,-1) or n ab = diag(-l, -1, +1).] 

It is well-known that the dynamics of class A Bianchi models derives from a Lagrangian which is obtained 
simply by substituting in the general action (|2.11[) the symmetry-reduced form of the metric, i.e. Eq. (|2.2p . This 
is also true when the metric is coupled to a homogeneous spinor. Here, we shall define the spatial homogeneity 
of a spinor VE' simply as meaning that the components of the spinor w.r.t. any frame h a :— h% which is 
invariant under the homogeneity group G, i.e. £% a h a = 0, depend only on timqj- 

We know on general grounds [29] that the lapse, N, and the shift vector, N a , will enter the final Hamiltonian 
action S — J{pdq — H(q,p) dt) as Lagrange multipliers of, respectively, the Hamiltonian constraint, TL, and the 
diffeomorphism (or momentum) constraints, Ha. Namely, the Hamiltonian is of the general form 

H = J{NU + N a U a )n (2.20) 

where /1 = t 1 A r 2 A r 3 is the spatial volume density [in the co- frame (dt, r 1 , r 2 , r 3 )], and where 

n = V9(2 4 Gg-r ) = V5(X- 4 <-r ), (2.21) 
n a = Vs(2 4 G2-r a VVs(2X-r a )- (2.22) 



Here, g denotes the determinant of the spatial metric g a b (w.r.t. the spatial-coframe r a , see Eq. (|2.2I0 . 
4 Gp = 4 i?Jg — \ A R5 c p the spacetime Einstein tensor, and Tg the matter stress-energy tensor. The factor 2 
multiplying 4 G a in the equations above represents (87rG) _1 in the units we use where 167rG = 1. 

Knowing in advance the structure (|2.20j) . one can simplify the computation of the Hamiltonian by working 
in the special quasi- Gaussian gauge where 

N = y/g , and N a = 0. (2.23) 

[This gauge was found useful in many previous cosmological studies, see e.g. [TJ [10l E].] In addition, we shall 
assume (for notational simplicity) that we consider the dynamics of a (comoving) piece of a homogeneous 
universe which has a unit (comoving) volume 1 = J ji. This allows one to identify the total Lagrangian with 
the Lagrangian density: 

S = fi dt L(q, q)= dt L(q, q) . 



4 Everywhere we use brackets [• • • ] around indices to denote antisymmetrization with weight one. 

5 In some cases, some spatial variation of ^, of the type £^ a 4* = i\ a NP, with real quantities A a subject to the integrability 
constraint \ a C a t c = 0, is compatible with the homogeneity of the geometry |16| . However, such a generalization is allowed neither 
in the case (we shall focus on) of a Majorana spinor, nor in the case of simple Lie algebras, such as type IX or type VIII, which are 
of most physical interest. 
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2.4 Gravity part of the Lagrangian 

The gravity part of the Lagrangian, which generically reads (modulo a total divergence) 

Leh = N ^g [g ac g bd K ab K cd - (g ab K ab f + R(g)] , (2.24) 
in terms of the spatial scalar curvature R(g) = 3 R(g), and of the second fundamental fornj^l, 

K ab := ^ (d t g ab -D a N b - D b N a ) (2.25) 

(D denoting the 3-dimensional covariant derivative), reads, when working in the gauge (|2.23|) 

L EH =T g (g,g)-V g (g). (2.26) 
Here T g denotes the "kinetic-energy" part of the gravity Lagrangian (g = dtg), 

T g {g, g) = \ g ac g bd g ab g cd - \ (g ab g ab ) 2 , (2.27) 

while V g denotes its "potential" part, 

V g (g) = -gR(g). (2.28) 

Using the decomposition (|2.6[) of g ab into the three diagonal variables (3i, /?2, f3 3 , and the three angles 
61,62,63 parametrizing the SO(3) [respectively 50(1,2)] matrix S a b in the type IX (resp. type VIII) case, we 
can express the gravity part (|2.26|) of the Lagrangian in terms of (3 a , 6 a and $ a ,6 a . 

The spatial scalar curvature of a general homogeneous metric g a b(t) T a (x) r b (x) is expressible in terms of 
the structure constants of Eq. (|2.3[) . namely (see, e.g., [3D]) 

R{g) = —-^C a bc C~hc — —C a bc C b ac — C a ac C a ac (2.29) 

where it is understood (for notational transparence) that a summation over, say, a and a denotes an appropriate 
contraction by means of g a b or its inverse g ab (e.g. A c Bc = g cc A C B C >). In the class A case the last term in 
Eq. (|2.29[) vanishes, while the other ones can be expressed in terms of the (symmetric) dual, n ad of C a bc (see 
Eq. (|2.19[) . This leads to the following simple expression for the "gravity potential" V g , Eq. (I2.28[) , 

V g (g) = n ab n^- 1 - (n a ^) 2 = g aa > 9 bb> n ab n a ' b ' - I (g ab n ab f . (2.30) 



Inserting the decomposition (|2.6j) into this result, then yields the expression of V g in terms of (3 a and the matrix 
S. As announced above, the fact that S a b was chosen to belong to the automorphism group A (leaving the 
structure constants C a bc invariant) implies that the right-hand-side (r.h.s.) of Eq. (|2.29[) is independent of S, 
and only depends on the diagonal variables (3 a . This is also true for the potential V g = — g R(g) if S belongs 
to the "special" subgroup of A having detS" =1. In the cases we consider here (types IX and VIII) this 
is automatically the case as A — SO(3) or 50(1,2). Finally, we conclude that V g (g) is given by the same 
expression that it would have if g ab had been assumed to be diagonal, namely 

V g {g) = V g ((3) = n\e~^ + n\e~^ + n\e~^ - ]- (me^ + n 2 e~ 2 ^ + n 3 e~ 2 ^) 2 

= I [ n 2 e -4ft + n 2 e -4& + n 2 e -4/33] 

- [ ni n 2 e- 2 ^ + ^ + n 2 n 3 e- 2 ^+^ + n 3 n ie - 2 ^ + ^} , (2.31) 
where n a denote the diagonal components of n ab = diag(ni, n2,n 3 ). 



3 Note that Ref. |29| defines K a f, with the opposite sign. 
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Turning now to the kinetic part T g , Eq. (|2.27p . we need to evaluate it in terms of the /3 — 5 parametrization 
(|2.6|) of gab- To be explicit we should, at this stage, choose a specific convention for the definition of the three 
(possibly generalized) Euler angles 6\, 62,63 parametrizing the (pseudo-)orthogonal matrix S. It is, however, 
better to introduce a notation for the "angular velocity", say w, of the matrix 5, before specifying its expression 
in terms of 61,62,63 and the 0's. We define w a g by writing 

S & c = J2 w& b st ' o ( 2 - 32 ) 

S 

or 5 = w 5, if we consider S = S a b and w = w a g as matrices. In this matrix notation, the decomposition (|2.6p 
reads g = S T e~ 2/3 5, where g denotes here the matrix g a b, (3 the diagonal matrix (3 a 8 a b, and the superscript T 
the transposition of a matrix. Differentiating the matrix g then yields 

g = S T {-2Pe- 2fi + e- 2fi w + w T e" 2 ' 3 ) S . (2.33) 

At this stage, the calculation depends on whether the matrix S is orthogonal (5*0(3); type IX) or pseudo- 
orthogonal (5*0(1, 2); type VIII). We shall henceforth focus on the type IX case, giving only some indications of 
the differences that arise in the type VIII case. In the type IX case we have (when using the usual basis where 
n ab _ jab^j gT g _ g gT — jj s0 that the "matrix angular velocity" w defined by Eq. (|2 . 32[) is antisymmetric in 
the usual sense: w T = —w. [In the type VIII case w would be antisymmetric in the Lorentzian sense, i.e. after 
considering :— rjaa> w a g where r/^ = diag(— 1, — 1, +1)-] Inserting this knowledge in Eq. (I2.33|) then yields 

an explicit expression for g a b of the form (denoting w ab = S b <= w a g = w a g) 

g cd = S & c S b d k-- b , (2.34) 

% = - 2 j3 a e- 2 ^ + (e- 2 ^ - e- 2/3b ) w &b . (2.35) 
Here, we can think of fc_^ as the components of the covariant tensor k a b '■= gab w.r.t. to the "rotated" co-frame 

T 7r :=5 a fc r b . (2.36) 

The co-frame (]2 .36[) is intermediate between the basic co-frame r a (w.r.t. which ds 2 = —N 2 dt 2 + g a bT a r b ), 
and the orthonormal frame 6 a (w.r.t. which ds 2 — rj & ^ 6 a 6" 3 ). Indeed, in the co-frame (dt,r a ) we have a 
non-Minkowskian, but diagonal form of the spacetime metric, namely: ds 2 = —N 2 dt 2 + ^2 e~ 2/3a (r a ) 2 , with 

N = Vfetg~; = e -(ft+fe+ft). 

As the gravitational kinetic-energy term (|2.27[) is manifestly invariant under any linear change of basis of 
the co- frame r a , it can be rewritten as 

T g = \ f c 9 Td Kb \l ~ \ it 1 \if . (2-37) 

where g ab are the components of the contravariant metric w.r.t. the intermediate frame f a , namely: g ab = 
e +2^a $ ab ^ This yields the explicit result 

T g =Tp0)+T w (w,P), (2.38) 
where ^ 

T p {P)=Y.fc- (E^) =-2(/3i /3 2 + /3 2 /3 3 + /3 3 , (2.39) 

a \ a / 

and (in the type IX case) 

Tl x (w,(3) = 2sinh 2 (/3 1 -/3 2 )( W i2 ) 2 +2sinh 2 (/3 2 -/3 3 )(u. 23 ) 2 

+ 2sinh 2 (/3 3 -/3i)(w 31 ) 2 . (2.40) 
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Summarizing so far: the gravity part of the Lagrangian has the form 

L EH (g,g) = T p 0)+T w (w,0) - V g {f3) , (2.41) 

with a p3- kinetic energy given by (|2.39p . a rotational kinetic energy linked to the dynamics of the Euler angles 
entering S(0 a ) given by (|2.40[) and a potential energy given by (|2.31|) (with n ab = S ab ). Note that the matrix 
S, i.e. the Euler angles 6*1,6*2, 6*3, do not explicitly enter the result (|2.41l) . As explained above, they do not 
enter the potential term V g because S was chosen to belong to the automorphism group of the Lie algebra ©. 
However, they do implicitly enter the rotational kinetic energy T w , as the rotational angular velocity w ab do 
depend both on 8 a and 9 a . For instance, if we define the Euler angles as in the standard references [3TJ I3"2"] . 
namely using the z — x — z convention, i.e. a matrix 

(cos ip smip 0\ (\ \ / cost/? smip 0\ 

-sin^ cosV' cos ^ sin0 J -sin</? cost/? o|, (2.42) 
1/ \0 -sinfl cos6y \ 1/ 

the angular velocities w ab (with w — S S -1 ) are given by 

w 12 = ipcosO + ip, 

w 23 = ip sin 9 sin ip + cos ip , 

w 31 = if sin 9 cos ip — 9 sin ip . (2.43) 

As several authors (see, e.g., |10j ) have previously remarked, the rotational part T w of the gravity Lagrangian is 
analogous to the kinetic energy of a rotating body ("asymmetric top"), namely T w = £l\ + I2 + ^3 ^3), 
where 0,\ = w 23 , f2 2 = w 31 , f2 3 = w 12 are the body-frame components of the angular velocity and J3 = 
4 sinh 2 (/3i — p}?), etc. the (body-frame) moments of inertia. Note, however, that, contrary to a usual rigid solid, 
the body-frame moments of inertia are time-dependent. Indeed, they depend on the /3's, which have a (coupled) 
dynamics determined by their kinetic energy Tp, their potential energy Vp, and their couplings to the other 
variables [and notably the Euler angles themselves, through the term T w (0 a , 9 a , p3 a )}- 



2.5 Spinor part of the Lagrangian 

Let us now derive the explicit form of the spinor part, Ld, Eq. (|2.13[) . of the Lagrangian. To do this, we need 
the explicit expression of the connection coefficients, say w^^. = ^> & p^ h * , of our specifically chosen vielbein 
9 a = h a ^dx^, defined in Eq. (|2.5[) in terms of the metric degrees of freedom /3 a , S(6 a ) (parametrizing g a b via 
Eq. (|2.6|) ). [In this subsection we take as above a vanishing shift vector iV a = 0, but allow for an arbitrary lapse 
N(t).} If we consider the "structure constants", say C a of the orthonormal co-frame 9 a , defined as 

M & = \C & (2.44) 
the (frame) connection coefficients w^^- can be expressed as 

where we denoted := r\ & & C a ^. 

As 9° — N(t) dt only involves t, we have d9° = 0, so that 

C% 7 = (and C m = Q). (2.46) 

On the other hand 9 a involves both t and the spatial coordinates (that are implicit in T a (x)). Using Eq. (|2 . 3f) . 
one finds that the Cartan differential of 9 a reads (no summation on a, but summation on b,c,d)) 

d9 & = d t (e-^ S & b )dtAr b + ^ e-^ S a d C d bc r b A r c . (2.47) 



Rewriting the r.h.s. in terms of 6° = Ndt and 8 a = e /3a S a b T b yields the structure constants C a and C a g g , 
namely (no summation on a, b, c) 

1 m" - 

C & 6i = ~C & M = "jy & 5 ab + e-"-+* (2.48) 

C & u = e-P' + ^'7f Vs , (2.49) 

where the C coefficients are the structure constants w.r.t. the intermediate frame T a = S a b T b , as defined in 
Eq. (12.91) above. Again, the choice of the matrix S as belonging to the automorphism group of the Lie algebra 
© implies that the C coefficients are simply equal to the original structure constants C a be- 

Inserting the results (|2.46|l . (|2.48l) . (|2.49p in Eq. (12.45[) (and considering the special case of type IX, i.e. 
w T = —w) yields the explicit expressions of the connection coefficients, namely 



1 r„ „ i 1 



otc = 2 ^ + °^ = N Ph 6bc + sinh(A - ^ ~N (2 ' 50) 

1 w ab 
= 2 [-°adb + Cbo & ] = - cosh(/3 a - ft) — , (2.51) 

uj hh& = \ (e-^+^+A C Q &c + e -ft+A>+/». C b ca - e-^'+P" C c ab ) . (2.52) 
In the last expressions, one has simply C a bc = £abc when using the usual basis r a for type IX, so that 

Lo &b£ = \ efc+A+/», (e~ 2 ^ + e" 2 ^ - e" 2 ^) s abc . (2.53) 

When inserting these results into the Dirac Lagrangian (|2 . 1 3[) . i.e. 

Ln=N^(y^d t y + *^L> &$ ^V-mmj , (2.54) 

there arise several types of terms: (a) a term involving <9 t (b) some terms involving dt f3; (c) terms involving 
the rotational velocities w ~ 9] (d) terms involving the purely spatial components u)~i - of the connection; and 
(e) the term involving the mass m of the spinor field. Let us first note that the lapse cancels out in all the 
terms involving one time derivative (dt dt (3 or w), while it contributes a factor N +1 in the "potential" terms 
(d) and (e). Let us first focus on the terms (a) and (b), i.e. those involving either dt "J/ or dt j3. They are easily 
found to be 



= 0**76^ + ^1^+..., (2.55) 

where we introduced the logarithmic derivative of gi = e~^^ 1+l32+ ^ 3 \ This shows (as had been used in many 
previous works), that the replacement of the original spinor variable \1/ by the rescaled spinor 

X ■= gi (2.56) 

disposes of the coupling to the $ a 's. This rescaling has also the effect of absorbing the prefactor ^fg in Eq. (|2.54[) 
into the various spinor bilinears: ^Jg^i{. . = x(. ■ -)x- Note, however, that the lapse prefactor N remains in 
factor of the bilinears of types (d) and (e), i.e. those involving no time derivatives. 
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Finally, after using the rescaling (12.56[) (and remembering the convention (|2.15l) ). we end up with a spinor 
part of the Lagrangian of the form 

L D = i X t X-cosh(/3 1 -/? 2 )w; i 2E i2 -cosh(/3 2 -/3 3 ) W 23 E 23 

- cosh(/3 3 - Pi)w %i £ 3i - U sgrav - U smass , (2.57) 
where we have introduced the short-hand notation 

S aS :=ix7 6 7 aS X=5X f 7 aS X (2.58) 

for the spinor bilinears that couple to the rotational velocities w ab (we used (|2.15p in the last equation). We 
are including a factor ^ in the definition (|2.58p so that the hermitian operator T, ab (in the quantum theory) 
measures the (second quantized) spin of the spinor field X (with eigenvalues ± ^ or 0; see below). 

In addition to the (body-frame) "spin-angular-velocity" coupling terms oc w ab E ab in the spinor Lagrangian, 
there are also (velocity-independent) "spinor potential terms'tj, that are naturally divided into two separate 
contributions: (i) the spinor potentials T4g rav coming from the coupling to the spatial connection coefficients 
w aSc' an d ( n ) the spinor potential U S mass coming from the Dirac mass term m^ty. The original expression 
(from (|2.54p ) of the gravitational-spinor potential is 

N - -s 

^grav^-^XTXs^X- ( 2 - 59 ) 

Using the gamma identity (where 7 a&c = 7 I a 7 b 7 c ]) 

and the fact that all the traces of w.g. vaniSn (because of the vanishing of C a a b), the potential (|2.59[) can be 
rewritten, using (|2.45l) . as 

V ssrav = -J2 jxr i£l &hd x- (2.61) 

a,b,c 

Using N = y/g — e -{0i+02+P3) anc [ (|2.49j) . this yields, in any (class A) Bianchi type (with n ab = diag(ni, n 2 , n 3 )) 

Kgrav = -\ me' 2 ^ xi^ X~\ n 2 e- 2 ^ X7 2M X - ~ n 3 e~ 2 ^ x^ i2 X , (2.62) 

where, for instance, the first term comes from the C 1 23 structure constant, and is associated with the corre- 
sponding (1; 2, 3) gravitational wall in /3-space, namely u>i ;2 3 = /3i — 02 ~ 03 + Pa — 2 /3i (when considering, 

a 

as we do here, the 3-dimensional case). We see (in agreement with Refs. that the corresponding spinor 

coupling involves 

Xi m x = iX i Toi m X- (2.63) 

Another peculiar feature of 3 dimensions is that all the different gravitational walls, wf. 23 , ^2-31 an d w 3-i2> 

involve the same spinor bilinear X7 123 X — X7 231 — X7 312 - This implies that, for instance in the type IX case, 
the spinorial-gravitational coupling explicitly reads 

K^av = ~\ (T** + e' 2 ^ + e~^) xi m X, (2-64) 
Finally, the last term in the spinor Lagrangian (12.571) reads 

V smass = +mN y fg V* = mNxx = me-^ + ^+^xx, (2-65) 



7 The classical analogs of these terms were discussed in Ref. 1171 . where x was treated as a c-number. 
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with 

XX = iX i 75 X- (2-66) 

Contrary to what happened for V s grav , this term does not correspond to one of the gravitational walls entering 
V g (j3). It would correspond to the wall form w\(j3) associated to a cosmological constant A. Indeed, in the 
N = y/g gauge, A generates a term oc AN y/g — Ae -2 "'*'^ with wa(/3) = J2 Pa- As usual (see, e.g., Ref. [21] ) 

a 

one expects the corresponding "spinor wall" to have a halved exponent, i.e. oc e~ WK ^\ which is the case of the 
spinor mass term (12.651) . 



3 Classical Hamiltonian formulation of the Einstein-Dirac Bianchi 
system. 

Having obtained a gauge-fixed Lagrangian formulation of the Einstein-Dirac system, let us now show how to 
pass to a Hamiltonian formalism. 



3.1 Prom the Lagrangian to the Hamiltonian in presence of derivative couplings 

A (well known) peculiarity of the Einstein-Dirac Lagrangian is the presence of derivative couplings between 
gravity and the spinor. In our homogeneous-cosmology context (and after the rescaling (|2.56[) of the spinor) 
these derivative couplings are the terms oc w ab E ab in (|2.57p , where w ab are linear in the time derivatives of the 
Euler angles (see Eq. (|2.43|l ). A well-known instance of such derivative couplings is given by the coupling of a 
charged particle to a magnetic field say (suppressing indices) 

L(q,q) = ^mq 2 + eqA-eV. (3.67) 

Let us recall the effect of the derivative coupling eqA, when going from the Lagrangian to the Hamiltonian 
formalism. First, it modifies the relation between the (canonical) momentum and the velocity, namely 

dL 

p= — = mq + eA. (3.68) 
oq 

Second, e A cancels out when computing the energy (because "a magnetic force does no work" ) 

E(q, q) := pq- L = (mq + e A) q - ^mq 2 - eAq + eV = ^mq 2 + eV . (3.69) 

However, third, the e A coupling reappears in the Hamiltonian because one must replace q by its expression in 
terms of p: 



H(q,p) = [E(q,q)] i(p) 



1 9 

- mf + e V 



{P - £A) +eV. (3.70) 



This mechanism is easily seen to hold for any derivative coupling which is linear in time derivatives. [In a more 
general case the mass m in (|3.67[l becomes some g-dependent quadratic form, and the m _1 factor in H(q,p) 
becomes the inverse quadratic form.] The crucial end result is that the Hamiltonian is numerically equal to 
the sum of the original kinetic energy terms (without the velocity coupling term eq A) and of the potential 
energy, but with the replacement of the concerned velocities by their expressions in terms of their (e A-shifted) 
canonical momenta. 



3.2 Hamiltonian formulation of the Einstein-Dirac system 

We can now write down the explicit Hamiltonian of the Einstein-Dirac system, for homogeneous configurations, 
when using the gauge-fixed vielbein (|2.4[) . (|2.5[) (with zero shift vector). More precisely, the Hamiltonian action 
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density has the form 



^Ham = ^2 1T l3 a Pa J r^2 Pe « ^ + * X f X 

a a 

-NU{p,K fi ,e,pe,x\x)-N a Ha{P,irf i ,6,pe,x\x), (3.71) 

where we have introduced the canonical momenta ^p a ,Pe a , respectively conjugated to the three /3's, and to the 
three Euler angles 8 a , and where N denotes the rescaled lapse N :— N/y/g. We shall work in a gauge where 

N = 1 and iV a = 0, which makes the Hamiltonian in the Hamiltonian action, Lnam = PQ — H, simply equal to 
H. Note that H = y/gH where H is the usual Arnowitt-Deser-Misner Hamiltonian density entering Eq. (|2.21l) 

(H is a spatial density of weight +1, while T-L has weight +2). Note also that we did not introduce any notation 
for the conjugate momentum of x as the Dirac action is first order, i.e. already in pq — H form (with p = i x 
when q = %). Using the results above (notably in the last Section) the value of T-L is the sum of kinetic-energy 
and potential-energy terms: 

H = Tp(7r j9 )+T IO (/3,7r IBS x,X t ) + W) 

+ V, grav (/3, X, X f ) + V s mass (/3, X , X f ) • (3.72) 

The various potential terms are the same as written above: V g (f3) is given by Eq. (I2.31[) . V SS i av (P, XiX ) by 
Eq. (12.62[) . and V s m ass(/?, X: X*) by Eq. (|2.65[) . [One should use x '■— * X^ 76 to replace x m terms of x* ■] 

The /3-kinetic-terrrJf| is originally given (in any dimension, and for any time gauge, i.e. any value of N := 

N/V9) by 



T P 0) = l,G ab p a p b = i 
N N 



(3.73) 



which defines the /3-space metric G a b- [Though we put the a index on (3 as a subscript, one should think of it 
as a contravariant index f3 a .] After the rescaling (|2.56|) of the spinor field, j3 appears only in the /3 kinetic term, 
so that the conjugate momenta 7r / 3 a to the /3 a 's are given by 

Trp. = tG ab p b . (3.74) 

This leads to the following value for the j3 kinetic energy expressed in terms of irp 

Tp(n e ) = ^G ab 7r^7r 0b , (3.75) 

where G ah is the inverse of G ab , i.e. (in our 3 + 1 dimensional case) 



a \ a 



G^tt^^^V^-- 1 . (3.76) 



Note that, if we were in a space dimension d ^ 3, the coefficient of the second term on the r.h.s. of (|3.76j) would 
be j^j instead of \. 

It remains to discuss the rotational kinetic term T w linked to the angular motion parametrized by the three 
Euler angles 8 a . To simplify this discussion, we shall henceforth come back to using the time gauge (|2.23l) . as 
we did in Section 2. [The formulas above for the (3 kinetic terms were written in a general time gauge as a 
reminder of the various occurrences of N and ^fg in the kinetic part of the action.] As the rotational velocities 
w ab ~ 9 C are linearly coupled to x ™ a the spinor bilinears S ah (see Eq. (|2.57|l ). we must apply the result of the 



8 When discussing the /3-kinetic-term in the Lagrangian (L = T — V), we do not take out a factor N , as we do in discussing the 
corresponding term in the Hamiltonian. 
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previous Section to this term. To ease the writing of the Hamiltonian version of the rotational kinetic term T w 
it is convenient to introduce the following momentum-like variables 

dL 

i.e. explicitly (in the type IX case) 

7T wab = 4 sinh 2 (/3 Q - t3 b ) w ab - cosh(/3 a - p b ) S afc , (3.78) 

as deduced from Eqs. (|2.40p and (I2.57[) . Here, and henceforth, we have simplified the notation by dropping the 
carets over the ("flat") indices a, b on w ab and E ab . [Note that the partial derivatives in Eq. (|3.77l) are done 
w.r.t. the three restricted independent components w 12 , w 23 , w 31 , of w; the other ones being defined in terms 
of these by, e.g., w 21 := — w 12 , etc.] 

Solving Eq. (|3.78l) for the w's as functions of the ir w 's, then yields the following Hamiltonian version (~ 
(p — eA) 2 1 (2m)) of the rotational kinetic term 

T w (p,n w ,x,X f ) = - . 7571 — +cosh(/3 1 -/3 2 )E 12 ] 2 

8smh [p\ - fa) 

I r , , r. ~ s ^23l2 



8sinh 2 (/3 2 -/3 3 ) 



[7r w 23 + cosh(/3 2 - /3 3 ) £ 



' V„3i + cosh(/3 3 - /3i) E 31 ] 2 . (3.79) 



!sinh 2 (/3 3 -/3i) 



If one wanted to express the Hamiltonian % completely in terms of canonically conjugated variables, one 
should replace the non-canonical, momentum-like variables tt w by the linear combinations of the Euler-angle 
conjugate momenta pg a (see Eq. (|3.7ip ) that they represent. As n w w — wdL/dw = OdL/dO = pg9, the 
transpose of the matrix A appearing on the r.h.s. of (|2.43l) yields the link pg = A t tt w , namely 

p v = cos 9 Tt w i2 + sin sin tp tt w 23 + sin cos ip Tt w 3i , 
pg = cos ip ir w 23 — sin ip n w 3i , 

Pip = n w i2 , (3.80) 

whose inverse reads 

Kw 12 = Pip 1 

. p v - COS , , 

77^23 = Sinw ; — \- COS ippg, 

smO 

p^-cosQp^, 

n w 3i = cosw— : — 1 smwpg. (3.81) 

smO 

The canonical Poisson brackets {9 a ,9b} = = {pe a ,Pe b } and {6 a ,pe b } = S a b are easily found to imply the 
following brackets for the non-canonical variables n w : 

{^12,^23} = —7^31, 

{7T TO 23 , 7T W 31 } = — 7T W 12, 

{77^31,^^,12} = — 7r w 23 . (3.82) 

Note that these are the opposite-sign brackets, compared to those of a usual angular momentum vector: {L x , L y } 
= +L Z , etc. Indeed, in the analogy of (|3.79[) (without the spin term) with the Hamiltonian of an asymmetric 
top, i.e. 

T 2 t2 t2 

H ^=k + k + k- (3 - 83) 
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the quantities Li, £2,-^3 (which are analogous to n w 2a,ir w 3i,w w i2) represent the (non-conserved) body- frame 
components of the angular momentum, which have opposite-type brackets, namely {^1,^2} = —L3, etc. This 
change of sign is linked to the fact that the body-frame Li,L2,Ls are related to the (normal-bracket, and 
conserved) space- frame L x ,L y ,L z by the time-dependent matrix S effecting the rotation between the space- 
frame and the body-frame. 

For most purposes, the expression (|3.79j) of the rotational energy in terms of the tt w 's, together with the 
knowledge of the Poisson brackets (|3.82[) between the tt w 's, suffices to discuss the Hamiltonian dynamics of the 
system. 

Summarizing: the total Hamiltonian of the gravity-spinor system is given by Eq. (|3.72l) . together with 
Eq. [Tp(wp)], Eq. [T W (P,1T W ,X)}, Eq. Eq. [Kgrav(/9,X)] and Eq. (£55) 

Ws ma,ss{Pi x)\- The Poisson brackets between the various variables is given by canonical pairs appearing in 
the pq terms in Eq. (|3.71[) . [See also (|3. 821) .] When considering "classical" spinor variables, one should use 
an odd, Grassmann- valued %, and define a graded (anti-) bracket under which \ is conjugated to ix ■ The 
quantum case will be discussed in detail below. 

3.3 Diffeomorphism constraints 

For simplifying the calculations, we have assumed above that we were working in a gauge where the shift vector 
iV a vanishes. However, as recalled in Eq. (|3.7ip . when we relax this assumption we know that N a will simply 
enter as a Lagrange multiplier of the diffeomorphism constraint H a ~ 0, where (see Eq. (|2 . 22[) 

H a = V? (2 4 R° a - T a °) = -2 D b -V9T°. (3.84) 

Here, the tensorial objects refer to the spatial metric g a b {e.g. D c g ao = 0), and all indices are projected on the 
co-frame r a , and its dual. In addition, 7rg in denotes the kinematical part of the conjugate momentum ir ab of the 
metric g a b- In the electromagnetic model discussed above, this would be the mq part, without the e A shift in 
Eq. (|3.68[) . In our case, n£ b n denotes the part of ir ab proportional to the second fundamental form K, i.e. (in 
covariant form) 7r^ n = ^fg(K a b — g a b K c c ), when using the + g a b convention (|2.25p for the definition of K ao . 
Finally, the matter term in (|3.84l) comes, in our case, from the stress-energy tensor of the spinor field. From 
Ref. [IS] , the explicit expression of the r.h.s. of (|3.84|) for a Bianchi spacetime reads 

H a = 2 7r£f n gbc C c a d + - S d c C c a d , (3.85) 

where S d c := h d a h ic S aS . 

Inserting the decomposition h a f, — e~~^ a S a b of the dreibein (with g ca - = 5~t h a c h b j), and transforming all 

tensor indices to the intermediate frame f a = S a b T b , except the indices on E ob which must remain "flat" (i.e. 
referring to 9 a — e"^ a T a ), one finds that (13.85[) becomes 

H^ := S a A H a = J2 ( 2 ^„ ^ 20v + \ ■ (3-86) 

p,q 

Note that, because of the covectorial nature of H a , the definition of Hm involves the inverse of the matrix 
S a b- The automorphism property of the matrix S guarantees that the r- frame structure constants C entering 
(I3.86|) are numerically equal to the original, r-frame, structure constants C . In addition, the rewriting of the 
ir ab g a b term in the Hamiltonian action as y]wf)„ $ a + ir w i2 w 12 + ir w 23 w 23 +ir w 3i w 31 , shows that the components 

a 

^kin = ^kin are linked to Trjffi via 

^=2(e-^-e- 2 ^)wll, (3-87) 
and similarly for the cyclic permutations 23 and 31. This implies that Eq. (|3.86j) reads, say for m = 3, 

% = - + cosh(^ - (3 2 ) Y} 2 . (3.88) 
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From Eq. f|3 . 78[) . we see that the last, spin contribution in (|3.88|) has the correct magnitude and sign to 
transform the kinematic contribution = 4sinh 2 (/3i — fc) w 12 into the full momentum- like variable ir w i2. 
Finally, this shows (in agreement with |18| ) that H-g = — tt w i2, i.e., more explicitly 

S a m Ha = — ir w pq , (3.89) 

where mpq is a cyclic permutation of 123. In other words, in the analogy where Tt w 23 — L\, etc. are the 
body-frame angular momenta, we have that "Hi = — L x , etc. are (minus), the space-frame angular momenta. 
As in the asymmetric top situation, the Ha, contrary to the ir w , are conserved because of the first-class bracket 
{HajH} — 0. Let us also note that the diffeomorphism constraints satisfy the Poisson-bracket algebra 

{H a ,n b }= -C c ab H c . (3.90) 

[Here, the minus sign comes from the minus sign in the relation "Hi = — L x .] The simple link (|3.89p between the 
diffeomorphism constraints and the (body-frame) angular momenta n w is not an accident but derives from basic 
symmetry properties. Indeed, while the three quantities H a generate (in phase space) adjoint transformations of 
the homogeneity group G |16j . the three angular momenta n w generate, by definition, the group of the matrices 
S(9i,02, 63) used in the parametrization (|2.5j> . (|2.6|) . However, for simple groups, such as type IX and type VIII, 
the latter group coincides with the adjoint representation of G. 

3.4 Explicit forms of some (class A) Bianchi Hamiltonians 

The full Hamiltonian action ()3.7ip implies two types of constraints: (i) the diffeomorphism constraints TA. a ~ 
linked to the arbitrariness of the shift vector N a , and (ii) the Hamiltonian constraint H linked to 
the arbitrariness of the rescaled lapse N = N/^fg. When considering general (class A) Bianchi models, the 
expression and the number of effective diffeomorphism constraints strongly depend on the structure constraints 
of the homogeneity group. Let us summarize the results for the various Bianchi types. 

In the Bianchi type I case, i.e. when C a be = 0, we sec on the general expression (|3.85p that Ha = 0, 
i.e. that there are no diffeomorphism constraints. In that case, one can still conventionally decide to use the 
decomposition (|2.6I) with a matrix 8(81,62,03) € £0(3). One then ends up with an Hamiltonian action of the 
form (|3.71[) . except for the last term N a H a which is absent. This type I action describes the coupled dynamics 
of the variables 02,03, 6\, 82,63, and x- The explicit form of the Hamiltonian is given (when N = N/ \/g) by 
the r.h.s. of p.72|) . where, however, several terms vanish because of the vanishing of the C"s. Specifically, the 
two potential terms in (|3.72[) linked to gravitational walls, namely V g (fi) and V sgTav ((3, x), are identically zero 
in type I. Finally, the type I Hamiltonian has the form 

U 1 =T {K )+T w {p,ir w , X ,X ] ) + V smass {f3,x,X ] ) (3-91) 

where Tp(iTp) is (universally) given by Eq. (|3.75[) . F smass (/3, x) by Eq. (|2.65l) . and where T w (f3, ir w ,x) is given by 
the full expression (|3.79j) . Note that this Hamiltonian (which must be submitted to the single constraint T-L ~ 0) 
describes the coupled dynamics of the variables /?i,/?2,/33; <p,6,ip, and x- It is a generalized asymmetric top 
dynamics (for the rotational degrees of freedom ip, 8, V>) , where the moments of inertia I3 (fi) — 4 sinh 2 (/?i — /? 2 ) , 
etc. have their own (coupled) dynamics, and which includes a coupling between the (bosonic) rotor and spinor 
degrees of freedom. As our main physical focus here concerns type IX, we shall not discuss in detail the type I 
dynamics. Let us only note that Ref. |23j has explicitly solved the equations of motion of type I dynamics, 
when using a different way of fixing the rotational state of the dreibein, and treating x as a classical quantity. 

The Bianchi type II case cannot be directly described by the expressions derived above from the SO(3) 
parametrization (|2.6j) because (contrary to the type I and type IX cases) its structure constants are not invariant 
under the full SO(3) group. Refs. [TBI H3 I2H] has indicated other useful choices of three-dimensional matrix 
groups leaving invariant the C"s. To give a different perspective on the type II case, and relate it to recent work 
on cosmological billiards, we treat the type II case in an Appendix by using an Iwasawa decomposition of the 
metric, i.e. an upper triangular matrix S in Eqs. (|2.5j) and (|2.6|) . 
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Leaving aside a discussion of the other degenerate class- A Bianchi types (i.e. VIo and VIIo), let us come 
back to the main focus of this work, namely the non-degenerate types IX and VIII. In both cases, the link 
(|3.89[) between the three diffeomorphism constraints, and the three body- frame angular momenta % w , holds. [As 
mentioned above, in type VIII the rotational angles and momenta refer to SO(l, 2) with metric diag(+l, —1, — 1), 
rather than to SO(3) with diag(+l, +1, +1).] We can therefore greatly simplify the dynamics by reducing the 
Hamiltonian by imposing the (first class) constraints T-L a — 0, i.e. tt w = 0, thereby eliminating both the pg 
terms in Eq. (|3.71[) and the N a *H a contribution. This leads to a reduced Hamiltonian action of the simpler 
form 

C = A, + iX f X - NH vm > IX W, tt-, X \ X) , (3.92) 

a 

where (Z = VIII, IX being a label for the Bianchi type) 

H z :=T p (n p ) +T^(p, X ) + V g z (p) + V s % av (p, X ) + V suiass (p, X ) ■ (3.93) 

Here: Tp(-Kp) and V smass (/3, x) have the universal structure given above; V z ((3) is given for any Bianchi type 
(including degenerate ones) by Eq. (|2.30[) ; V s z grav is given for any Bianchi type by Eq. (|2.61[) , with C a g e given 
by Eq. (|2.49j) (with C = C when S leaves the C's invariant, as is the case for Z = VIII and IX); and, finally, 
the rotational energy term takes, in type IX, the simplified form obtained by replacing n w — ¥ in (|3.79|) 

T<° )IX (/3,X) = icoth^x-AKS^ + lcoti^-AXS^) 2 

+ icoth 2 (/3 3 -/?i)(S 31 ) 2 . (3.94) 

The type VIII case is obtained by particularizing to the case (ni, n 2 , n 3 ) = (+1, +1, —1) the result for a general 
n ab = diag(ni, ri2, n.3), which reads 

1 /„, P Pa-Pi 4. noP ft-&\ 2 

^ d - 5 tz-» ) <e > + cycuc ■ <3 ' 95) 

The Hamiltonian (I3.93p ~ p.95p describes the coupled dynamics of the diagonal degrees of freedom of the 
metric coupled to a spinor. In the next Section we shall discuss its quantization. 



4 Quantum formulation of the coupled spinor-Bianchi-IX system 
4.1 Dependence on the Euler angles 

We shall denote by $ the abstract quantum state of the system. Following Dirac, we interpret the various 
classical constraints, T~L rs 0, H a ~ as constraints on $ of the type 

H$ = Q, (4.1) 

£a$ = 0, (4.2) 

where % and T-L a are suitably defined operatorial versions of % and "H a . We shall work in the gravity configuration 
space Pi,^2,^3] Oi,02,03, together with a suitable description of the spinor degrees of freedom, labelled, say, by 
a (see below). This leads to a wave function of the universe described by < I > (/3, 0, a). Let us start by remarking 
that the classical link (|3.89l) naturally suggests an ordering such that Eq. (|4.2p becomes 

n wab $(p,0,a) =0. (4.3) 

The classical angular momenta tt w are linear combinations of the canonical angular momenta pg a , see 
Eq. (|3.8ip . The canonical quantization p$ a = —id/d8 a (with H = 1), together with the natural ordering 
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of the tt w as differential operators on the 50(3) spac^l naturally leads to concluding that the quantum con- 
straint (|4.3j) means that #2, #3) does not depend on the Euler angles: de a & = 0. More directly, we can 
say that Eq. (|4.3j) means that $ is a zero-angular-momentum state on SO(3); i.e. an "s-wave". As $ does not 
depend on the 0's, we can simply ignore them in the following. Note that this is consistent with having started 
the quantization procedure at the level of the final, reduced type IX Hamiltonian, namely Eq. (|3.93[) . which 
does not contain the (9's. From this point of view, the quantum description is cleaner than the classical one. 
Indeed, in the classical description the Euler angles satisfy first-order differential equations obtained by setting 
to zero the l.h.s. of Eq. (|3.78l) . This yields a non-trivial dynamics for the 0's sourced by the spinor bilinears 
oc and influenced by the time-dependence of the /9's. In the quantum description, the 9's do not appear 
at all, which is nicely consistent with the fact that they have the character of gauge parameters. 



4.2 Quantum description of the spinor degrees of freedom 

After the rescaling (|2.56[) of the spinor, the spinor kinetic term is simply 

4 

Lf» = i X h = iJ2xl(t)x a (t) (4.4) 

Q = l 

where we made explicit the Dirac-spinor index a = 1,2,3,4, which was kept implicit up to now. In the case 
where \ lS a generic Dirac spinor, with 4 independent complex components, the quantization of (|4.4I) is done 
by the standard anticommutator result (with h = 1) 

[Xa,Xp]+ = = [xt,Xp]+; [Xa,Xp]+ = <W , (4.5) 

where [.A, !?]+ = AB + BA denotes an anticommutator, and where all the x' s are taken at the same instant 
t. In other words, each one of the four x» s can be viewed (at each instant t) as an independent (complex) 
fermionic destruction operator, with x a being the corresponding fermionic creation operator. To each pair 
(Xa, Xa) then corresponds a two-dimensional (complex) Hilbert space, so that to the 4 mutually anticommuting 
pairs (xa, Xa)i a = ^1 ■ ■ ■ >^ correspond a fermionic Hilbert space of dimension 2 4 = 16. 

We can, however, simplify the problem by demanding, from the start, that x be a Majorana spinor, i.e. 
be restricted to contain only 4 independent real components. More precisely, let us use the following explicit 
Majorana representation of the 4 Dirac gamma matrices 7" entering the spinor action: 



i (72 \ T _ ( o~i \ 2 _ f ~ °3 



[10-2 ) ' T I. g x ) ■' \ 

where o~\ = ( ^ ^ , 02 = f ^ , 03 = ^lj &Te ^ G s ^ anc ^ ar< ^ P aun matrices. The four gamma matrices 

(|4.6I) are all real, with the 3 spatial 7's being symmetric (and therefore hermitian), while 7 is anti-symmetric 
(i.e. anti- hermitian) . 

When using a real representation such as (I4.6[) . the reality condition for a Majorana spinor is simply x* = X> 
i.e. the condition that each component Xa be real (i.e. hermitian, as an operator). In that case, the quantization 
of the standard spinor kinetic term (|4.4[) . i.e. iJ^XaXa, involves an extra factor i, namely 

a 

[X<x,X0]+ =XaXP+XpXa = ^<W- (4.7) 



9 Note that the 8's parametrize SO(3) matrices, rather than the SU(2) group G which constitutes the cosmological space. The 
"factor 2" difference between SO (3) and SU(2) might have been important, if we had needed to consider situations where tt w 4> ^ 0. 
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The need for the factor | in the quantization condition (|4.7[) can be viewed in several ways. The most direct 
way is that the normalization of the basic field commutators should be chosen so that the universal Heisenberg- 
representation equation of motion for operators (with h = 1) 

ij t Q={Q,H] (4.8) 

(where [A, B] = AB — BA denotes a commutator) hold true, and be consistent with the usual Euler-Lagrange 
equations of motion. As the Euler-Lagrange variation of the kinetic term 8 J dtixx (for a real, anticommuting 
x) involves a factor 2 (after integrating by parts), one needs the factor | in (|4.7|) . 

An indirect check consists of decomposing a general, complex Dirac spinor, with kinetic term (|4.4[) . into real 
(i.e. Majorana) and imaginary parts, say x — XR. + i Xi, X 1 — XR ~ * Xii where XR an d Xi are both Majorana. 
It is then easily checked that the standard anticommutators (|4.5I) imply that xr an( i Xi both satisfy (|4.7p . and 
mutually anti-commute. In addition, it is easily seen that if we insert the decomposition x = XR + i Xi m the full 
Dirac action (|2.57l) . it simply decomposes as the sum of two decoupled Majorana actions, one involving xr. an d 
one involving xi- [This is true because all the spinor bilinears, including the mass term oc X^ToX, involve an 
antisymmetric Clifford- algebra matrix, such as 7 1 j 2 , 7g 7 123 or jq, sandwiched between x* an d X-] This shows 
that considering a complex Dirac spinor source is equivalent to coupling gravity to the sum of two independent 
Majorana spinor sources. At the Hilbert-space level, the Dirac fermionic state space is just the tensor product of 
two independent Majorana fermionic state spaces. For simplicity, we shall consider in the following the simpler, 
single Majorana case, i.e. x* — X- Note, in particular, that in that case the fermionic Hilbert space is simply of 
dimension 2 2 = 4. Indeed, this fermionic Majorana space is the (irreducible) representation space of the algebra 
(|4.7p . which is simply a Clifford algebra on the 4-dimensional Euclidean space. Actually, if we introduce the 
notation 

T a :=2 Xa , a = l,...,4 (4.9) 

we see that the four operators (in Hilbert space) T a , a — 1, ... ,4, satisfy the standard 0(4) Clifford algebra 
relations 

r a r> + i>r Q = 26 a0 . (4.10) 

In particular, each r Q has a unit square. In addition, note that the product 

T5 := rir 2 r 3 r 4 = 16x1x2x3x4, (4.11) 

defines an operator which anticommutes with each Xa, a nd which has also a unit square: r 2 , = +1. [This 
contrasts with the 0(3, 1) case where 75 = ToTf TjTi squares to minus 1, because of the time-like character of 
75-] 

Summarizing: when x is taken to be a Majorana spinor, the second-quantization of x, Eq. (|4.7p . leads to a 4- 
dimensional space of fermionic degrees of freedom. In other words, the wavefunction of the universe, $((3 a , 8a, c), 
is labelled, besides the 6 continuous bosonic labels (/?i , (3 2 , /?3 ; Oi, 62,83) of the gravitational minisuperspace, by 
a discrete index a, which takes 4 values, and which is algebraically equivalent to the spinor index of the Dirac 
representation of the 0(4) Clifford algebra (|4.10j) . When x is taken to be a complex, Dirac spinor, the second 
quantization of x, Eq. (|4.5p . leads to a 16-dimensional space of fermionic degrees of freedom, which is the tensor 
product of two independent Majorana-spinor spaces. Note that if one works in a Heisenberg picture, where the 
x(i)' s evolve in time according to the equations of motion (|4.8p . the corresponding Clifford algebra generators 
(|4.9p must also be viewed as time-evolving objects, rather than fixed numerical gamma matrices. 

4.3 Explicit form of the Einstein-Majorana Bianchi-IX quantum Hamiltonian 

After having taken into account the diffeomorphism constraints (|4.2p in the form (|4.3p , so that the wave function 
$(/3, cr) depends only on the three diagonal metric variables /3i,/3 2 ,/33 and the discrete spin-space label tr, it 
is natural, in the gauge N = N/^/g = 1, to quantize the gravitational degrees of freedom by replacing the 
classical (3- momenta up [which enter the Hamiltonian (13.931) only through Xjg (71-/3), Eq- (|3-75p ] by the differential 
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operators iip a :— —id/d/3 a . This leads to interpreting the Hamiltonian constraint as the following Klein- 
Gordon-like Wheeler-DeWitt equation 

H Ix ^(l3,a) = 0. (4.12) 

with ^ 

H 1X := T p (n ) + f^(0, X ) + V^ x ((3) + \>^ rav (/3, x) + V smass ((3, X ) + c , (4.13) 
W/O = - \ G ab 3 Pa d Pb = -\ (dl + 81 + d%) + 1 (5 A + 9ft + 9ft) 2 , (4.14) 

where Vg(/3) is the usual type- IX potential (|2.3ip . and where the quantum versions of all the ft — x coupling 

terms Tx°\/3, X )i V^. av (/3, x)> Kmass(/3)X)) are simply obtained from their classical expressions by interpreting 
the various spinor bilinears they contain as operators in the fermionic space defined, say in the Majorana case, 
by Eq. ([377]) . This yields 

Tl°Hf3, X ) = gCoth 2 (A-&)(£ 12 ) 2 + ^coth 2 (^-/3 3 )(S 23 ) 2 

+ l co th 2 (/3 3 -/3 1 )(S 31 ) 2 , (4.15) 

where 

£ 12 :=^xV ? X, S 23 :=^xV 3 X, S 31 := | xV T X , (41g) 

are operators acting in the 4-dimensional fermionic space. Similarly, denoting by "+ cyclic" the addition of two 
other cyclically permuted terms, 



V^rav(Ax) = - \ ^ 2ft (u f 75 7 T5S X)+ cyclic 

= -Ue-^+e-^+e-^xhol^x), (4.17) 



and, 

V smass (/3,x) = me-^+^ + ^(i X ' f 7oX)- (4.18) 

Finally, we have allowed for the addition of an "ordering constant" c in Eq. (|4. 13[) . There are several motivations 
for allowing for such an additional constant. First, if we were working in a different gauge, e.g. N = 1, i.e. 
N = g- 1 / 2 = e (0i+#>+#O instead of N = 1, we would have had to quantize H = \ e^ 1+ ^ 2+M G ab %p a np b 
instead of jG ab wp a 7r^ b . More generally the quantization of % = | e 2 ™^' G ab 7Tft 7Tft , where a(/3) = a a f3 a is 
some linear form in the /3's, leads (after reabsorbing e"'^ in the definition of the wavefunction) to an ordering 
ambiguity which is equivalent to adding a constant c oc G ab a a ab in Eq. (|4.13[) . Second, the contribution (|4.15l) 
being quartic in the non-commuting X 's has also ordering ambiguities (though it is natural to define it as we 
did, i.e. by taking the squares of the spinor bilinears (|4.16p . which have no ordering ambiguities because 7 12 , 
etc., are antisymmetric matrices). Finally, there is a basic ambiguity in the gravity-spinor Lagrangian due to 
the possible use of first-order (Palatini) versus second-order formulations. It has been shown long ago [33 that 
this leads to an ambiguity in the action density proportional to 

A£ = ^V TI 3(?7 3? ^)(?7 3? ^) 

= N[-(i X t x f - (j X t 7 2 3 x) 2 _ ( - x f 7 5T x) 2 

+ (*X t 7o7 T53 x) 2 ]. (4.19) 

We will come back below to the latter (/3-independent when N = 1, but x-dependent) ambiguity. Finally, 
note that one could also consider adding to the quantum Dirac action other Lorentz-invariant contributions 
quadratic in spinor bilinears, A£ ~ (ipj A ip) 2 , because no symmetry prevents the appearance of such quartic 
contributions. 
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4.4 Properties of the spinor bilinears 

For the time being, let us note that, in the representation (|4.6I) . one can compute explicit expressions for the 
Majorana-spinor bilinears entering the Hamiltonian. Let the 4 components of the Majorana spinor \ m the 
representation (|4.6[) be denoted as Xi?X2,X3?X4> an d let \ T denote the transpose of the column vector \a, i.e- 
the horizontal row \ T — (xii X2, X3i Xi)- We find 



t 1 2 

X 7 7 X 

T 2 S 

x 7 7 x 

X 7 7 X 



Xl X2 - X2 Xl + X3 X4 ~ X4 X3 , 
X2 X3 - X3 X2 + Xl X4 - X4 Xl , 
X3 Xl - Xl X3 + X2 X4 - X4 X2 , 



for the bilinears entering (without the i/2 prefactor), while 



X 70 71 72 73 X = Xl X2 - X2 Xl - X3 X4 + X4 X3 • 



and 



(4.20) 

(4.21) 
(4.22) 



X 7o X = X2 X3 - X3 X2 - Xi X4 + X4 Xi , 

are the bilinears entering V ssrav and V smass . 

In view of (|4.20|) it would seem that Tw\ which contains the squares of these bilinears, is a hopelessly 
complicated expression. However, things simplify very much if we use the Clifford-algebra properties of the 
X's. More precisely, using the notation (|4.9[) , and simple properties such as r 2 = 1, (ril^) 2 = —1, and 
Ti T2 T$ — —T3 T4 (with the definition (|4. 1 1[) of T$), one finds the following simplified expressions for the spinor 
bilinears 



T 1 ? 

X 7 X 


= rj 


r 2 


XV 3 X 


= r 2 


r 3 


X T 7 31 X 


= r 3 


ri 


X T 70123 X 


= r x 


r 2 


X T 7oX 


= r 2 


r 3 



1 


-r 5 




2 


1 


-r 5 




2 


1 


-r 5 




2 


1 


+ r 5 




2 


1 


+ r 5 



(4.23) 



Note that these expressions contain the projectors 



i±r 5 



(4.24) 



which satisfy the usual properties of Dirac's helicity projectors (1 ± i7s)/2 [with 75 = 70717273 so that 
(i 7s) 2 = +1], namely: P+ + P_ = 1, P 2 = P + , P 2 _=P_, P+ P_ = 0. 

Let us also note that the above bilinears satisfy simple commutation relations among themselves. First, as 
r a fc '■— r[ a r h ] (i.e. Tn — 0, Ti 2 = V\ T%, etc.) commutes with and P+ P_ = 0, we see that the first three 
bilinears in (|4.23[) commute with the last two. Second, one easily checks that (together with its cyclic analogs) 



^12 ^23 ' 
T J_ T 7_ 

X 2 X 1 X 2 X 



^13 
T ]_ 

X 2 X ■ 



(4.25) 



i. e. adding the factor i transforming the anti-hermitian \ % T '-f ab ^ into the hermitian operators S afc of Eq. (|4.16l) 



^12 ^23 



31 



(4.26) 
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together with its cyclic kins. In other words, the second-quantized operators E 12 , E 23 , E 31 satisfy the (reverse- 
sign) commutation law of body-frame angular momenta (the classical Poisson bracket {Li, L2} = —£3 mentioned 
above becoming [L\, L2] = —1L3 at the quantum level). 

Let us note that, more generally, one can deduce directly from the basic anticommutation relations (14.71) 
that the commutators of two bilinears x T Ax = Xa A aj s xp and X 7 B X = Xa B a p xp (where A and B are some, 
say antisymmetric, ^-independent matrices) is equal to 

[X T A X , X T B X }=X T [AB}x, (4.27) 

where [A, E\ a p = A aa B a p — B aa A a p is the usual matrix commutator. 

In addition the squares of the bilinears (14.23)) simplify. Notably, (x T 7 12 x) 2 — (I\ T 2 ) 2 Pi = -P_, i.e. 

(e 12 ) 2 = (E 23 ) 2 = (E 31 ) 2 = \p-- (4-28) 
4.5 "Squared-mass" in the Wheeler-DeWitt equation 

The main result of Ref. |23j was the finding that a classical spinor source (i.e. treating the spinor bilinears 

X T Ax as c- numbers) modifies the usual quadratic relation ^Zv\ — (12 Pa ) — satisfied by Kasner exponents 

a \ a J 

(in the Bianchi type I case, i.e. far from any potential wall) into the relation 10 ! 

(4.29) 



E^ 2 - (x>«) =-\ t( s12 ) 2 + ( s23 ) 2 + ( s31 ) 2 ] 

a \ a / 



where the E af> are the spinor bilinears defined above (treated as c- numbers). In the language of cosmological 
billiards (and using the time gauge N = Jg = T, i.e. a coordinate time t = — InT), the result (|4.29l) means 

that the Lorentzian squared velocity in /3-space /3 2 := G a bPaPb = J^Pa ~~ ( J^Pa ) is negative, p 1 < 0, i.e. 

a \ a / 

time-like, rather than light-like (p 2 = 0) as usual. Belinsky and Khalatnikov [23] argued that this implied 
that a classical spinor would (like a scalar field does) ultimately quench the chaotic, oscillatory regime near the 
singularity, and ultimately turn it into a monotonic Kasner-like, power-law behaviour (because the /3-particle 
can ultimately move on a time- like geodesic which no longer hits any gravitation wall). 

The result (|4.29[) was understood via an Hamiltonian approach, and within a more general Kac-Moody 
coset model approach, by de Buyl, Henneaux and Paulot [23] ■ These authors found that, far from the walls 
associated to a general coset model, the mass-shell condition for the /3-momenta 7r^ a cx G a b 0f, was modified, 

t 2 = (where 7T 2 , 



by the coupling to a classical spinor, from the usual light-like condition 7rS = (where 7r? := G ab irp a tt^) to a 



massive-shell condition 
where the general expression for the (/3-space) squared-mass ^ 2 was found to be 



M 2 = , (4.30) 



2 1 V.. ; , t I- > ,'-2 

coset 2 / j 



i^Jlxf- (4.31) 



Here, x is the coset Dirac field [which corresponds to the rescaled Einstein-Dirac field (|2.56p ]. and a labels 
all the (positive) Kac-Moody roots (including their degeneracy) that enter the bosonic coset model G/K. 
Each root a is a linear form in the /3's, and corresponds (in Iwasawa gauge) to a bosonic wall potential 
proportional to exp(— 2a((3)). The object J* in (|4.31|) is an anti-hermitian generator which represents (in 
the spinor representation s) the "rotation" E a — E~ a of the maximally compact subgroup K of G associated 



°We normalize the type I spatial metric so that yj~g = T, where T is the proper time. 
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to the root a. The extra factor i in (14.31| turns J* into a hermitian operator, so that i~xJ Ja X is rea l an d A* 2 is 
given (for a Kac-Moody coset G/K) by an infinite sum of positive quantities. 

In the language of the general Kac-Moody result (|4.30[) , (|4.31[) (which uses a normalization where 7r^ a = 
G a b$b, without the factor 2 entering our Eq. (|3.74[) above) the original Belinsky-Khalatnikov result (|4.29|) 
corresponds to the squared-mass 

/4 K = i[(S 12 ) 2 + (S 23 ) 2 + (S 31 ) 2 ] . (4.32) 

This squared-mass is of the general form (14.311) ( see below) but includes only three roots, namely the roots 
corresponding to the so-called "symmetry walls"! ai2(/3) = fi 2 — /3 1 , 023 (/3) = /3 3 — /3 2 and «i3(/3) = f3 3 — (3 1 . 
One might wonder whether the presence of only those three symmetry roots, and, in particular, the absence of 
any contribution coming from the gravitational- wall roots af. 23 = 2/3±, etc. (which are crucial to generate 
the BKL chaos) is due to the approximate nature of the treatment of Ref. [23] , based on a Bianchi type-I model, 
i.e. a model which contains all the symmetry walls, but no gravitational walls. Separately from this issue, we 
wish here to use our results above to clarify the effect of considering, as one should, the spin source as being 
quantum, rather than classical, as was assumed both in [33] and in [2~I] . 

One can define the squared-mass term in the Wheeler-DeWitt equation as the term remaining with tt^ 
in a BKL-type limit j3 — > +00 in which one stays far away from all exponential walls. We have already 
mentioned that Vj^. av , Eq. (|4.17j) . contains (half) gravitational walls, while V s mass, Eq. (|4.18j) . contains a (half) 
cosmological-constant-related exponential wall. All these terms, as well as the usual bosonic type-IX potential 
V lyi ((3), Eq. (|2.3ip . will tend to zero in this far-wall BKL limit. It remains to consider the contribution 
Eq. (|4.15p . Using the hyperbolic-trigonometry identity coth 2 x = 1 + 1/ sinh 2 (x), we see that we can split 
as 

^ 0) = 7^l + ^ CCntrif (^x), (4.33) 



4 ' 

where 

[(£ 12 ) 2 + (£ 23 ) 2 + (£ 31 ) 2 , (4.34) 



~2 1 



and 



o smh (/3i - f3 2 ) 

One recognizes in (|4.35l) the sum of three sinh _2 -type centrifugal-wall potentials, which is the form taken by 
symmetry walls when one uses a Gauss decomposition of the off-diagonal components of the metric rather than 
an Iwasawa decomposition (see [TU1[TT]). Far from the symmetry (or centrifugal) walls, i.e. when |/3 2 — /3 X | S> 1, 
etc., v ccntTl1 becomes a sum of exponentially small terms oc exp(— 2 \(3 2 — /3 1 1 ) + cyclic. This leaves as only terms 
contributing to the dynamics in the far- wall limit 

H IX = - [tt 2 + ^| + 4c] + (exponentially small terms), (4.36) 

where 7? 2 = —G ab dp a dp b is the Klein-Gordon operator associated to the /3-space Lorentzian metric. 

This result shows that, in Bianchi type-IX (and also, as one easily sees, in type VIII), the coupling of gravity 
to a quantum spinor generates a g-number squared mass in the Klein-Gordon-like Wheeler-DeWitt equation 
given by 

V 2 = Ms + 4c , 

where /If, is defined in Eq. (|4.34j) . and where c is the ordering ambiguity that we allowed for in Eq. (|4.13j) . 

Note, first, that if we take a formal classical limit of /If;, Eq. (|4.34[) . we recover the Belinsky-Khalatnikov 
result //| K , Eq. (|4.32[) . without any extra contribution coming from the gravitational walls (which were explicitly 
taken into account in our calculation, contrary to the type-I computation of Ref. [23j). When comparing this 



^When considering potential wall forms w^(/3) from a coset perspective we denote them as a^(fi). 
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result with the coset-model result /i 2 oset , Eq. (|4.31[) . we see that there is a genuine difference. Even if one were 
considering a truncated coset model involving only the couplings to the walls entering the type-IX gravity model, 
the sum in Eq. (|4.31|) should include as roots both the three symmetry roots a.12 (/?) = /fe — «23 {0) = ^3 — /?2 
and ai3(/3) = ,83 — /3i, and the three gravitational roots af 23 (/3) = Pi — P2 — 03 + S a Pa = 2 /3i, a^siiP) =2/32 
and afi 2 (/3) = 2/3 3 . The corresponding K(AE 3 ) generators, J a — E a — E- a , within the spinor representation 
of K(AE 3 ), are known to be El E3 ESI 122] , 
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7 s 



2 7 
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(4.37) 



for the three symmetry roots, and 



J «? 23 " 2 70 7 



7 s ~, „,231 

Q 231 _ 2 ' 



J a 9 = n 70 7 

"312 2 
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(4.38) 



for the three gravitational walls. Inserting the expressions (|4.37l) into the general coset result (I4.31|) yields, 
in view of the definition (|2.58l) of precisely the result (|4.34[) (including the correct normalization factor). 
However, a coset model involving the couplings to the walls entering the type-IX gravity model would yield a 
(quantum) squared-mass of the form 



~2(IX) 
P-cosct 



#1 



2 V 2 



,123 



(4.39) 



where the factor 3 comes from the fact that the three gravitational-root generators (I4.38P happen to be equal 
among themselves. 

The reason why the coset-results (|4.31[) or (|4.39p (for the type-IX truncation) contain more contributions 
than the original type-IX gravity model, Eq. (|4.34[) . is easily understood from the corresponding derivations of 
these results in [24] and in Section [2] above. Indeed, one sees that each individual contribution 5 (« x'-'aX) 2 
to /1 2 comes for the presence of a corresponding time-derivative coupling ~ eqA, see Eq. (|3.67[) . between 
gravity and the spinor. For instance, it is because in Eq. (|2.57l) there were (say in the far-wall limit) three 
Lagrangian-coupling terms, 



cosh(/3! - w 12 S 12 + cyclic : 



1 



e a 12 (/3) w 12 S 12 



cyclic , 



(4.40) 



between the rotational velocities w 12 = (ficosO + ip, etc. of the off-diagonal metric variables (parametrized by 
Eulcr angles), that one ended up with /if; given by the sum of three symmetry wall contributions ai 2 , 0^3, 013. 
By contrast, the couplings corresponding to the three gravitational walls af 2 3, ^231 > a 3i2' between gravity and 
the spinor, predicted by the Einstcin-Dirac action, were contained in V^grav, Eq. (|2.64[) . i.e. were of the type 



s grav 
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X f 7o7 123 X 



cyclic . 



(4.41) 



This is very analogous to the coupling (|4.40|) . with, however, the crucial difference that the rotational 
velocities w 12 , etc., appearing in (|4.40p are replaced in (|4.4ip by the structure constants C 1 2 3, etc. In the 
gravity picture, the structure constants are spatial derivatives of the metric, and therefore do not contribute a 
squared-mass term via the 



(p-eAy 
2m 



1 e- 

2 ~ 



A 2 



(4.42) 



mechanism (with eAoc^ X Ja X) explained above, which was linked to a time- derivative coupling AC = eqA. 
By contrast, the coset action is related to the gravity action by a dualization of some of the gravity degrees of 
freedom. In particular, the structure constants C a be become replaced by the conjugate momenta, i.e. essentially 
by the time-derivative, of some dual coset field, namely C ab := \e hcd C a d oc n afc ~ <p a t for the AE 3 case linked 
to the (3 + l)-dimensional Bianchi-IX dynamics (see Section 8 of [S])- This explains why the gravitational roots 
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(and actually all the coset roots) give a contribution to /j? within the coset model (for the spin 1/2 case), but 
not within the original gravity model. 

Let us now turn to the quantum nature of the Bianchi-IX squared-mass term (|4.34[) . Though each term 
(E° & ) 2 = (l x^ ab x) is quartic in the quantized spinor field x, we have shown above that the Clifford- algebra 
nature of the quantization conditions for a Majorana spinor led to simple final results. More precisely, we see 
from (|4.28[) that each separate (E af> ) 2 term yields the same result, so that the total quantum fi 2 reads 

where we recall that T$ denotes the involutive (r 2 = 1) operator f|4. 1 1|) . The simple result (|4.43[) shows that, 
in the quantum theory, the squared-mass term [in absence of any additional ordering-constant contribution 4c, 
as in Eq. (|4.36[) ] is an operator which has two possible eigenvalues: the eigenvalue ^ = in the 2-dimensional 
space of "helicity" T 5 = +1, or the eigenvalue = 3/8 in the complementary (orthogonal) 2-dimensional 
space of "helicity" T$ = — 1. [Note that ^ vanishes in half of the total (4-dimensional) Majorana-spinor 
Hilbert space.] In terms of (2j + l)-dimensional irreducible representations Dj of SU(2), we can view the 
three components E 12 , E 23 , E 31 of a single Majorana field as describing the sum Dq + Dq + Di. Let us also 

note that our results above, also gives the eigenvalues of E 12 ,E 23 ,E 31 (and then of /i 2 ) when considering a 
complex Dirac spinor. Indeed, this case is obtained by considering x = XR + * Xi as t ne complex combination 
of two, independent (i.e. anticommuting) , Majorana spinors, so that each final E 12 etc., is given by a sum, 
Ejf + E} 2 , of two commuting spin operators. The corresponding representation space is therefore the tensor 
product (D a + Dq + Di) x (Dq + Dq + Di), i.e. the sum 5D + 4Di + D\. This shows that the possible 
eigenvalues of jl 2 are 0, 3/8 and 1. Moreover, the eigenvalue /I 2 = is obtained (when assuming c = 0) in 
a 5-dimensional subspace of the 16-dimensional total Hilbert space. Note that we can extend this result by 
considering the case where gravity couples to a sum of N independent (i.e. anticommuting) Majorana spinor 
fields. In that case each spin operator E 12 etc. will be the sum of N independent (commuting) contributions 
belonging to a Dq + Dq + Di representation space. This implies that the eigenvalues of jug will range between 
and N(N + 2)/8. Decomposing the Hilbert space obtained from N tensorial products of (Do + Dq +Di) into 
a direct sum of irreducible spin s subspaces D s , one finds that the number A(s,N) of irreducible spaces D s 
appearing in the 4^ dimensional product space is given by : 

A(SjjV) _(2 S + D (2N + 2)! 



N + l (TV + 2s + 2)! (TV -2s)!' 

For large N and fixed s this behaves as : A(s, N) ~ 4 JV+1 (1 + 2 s)/(y/Tr N 3 / 2 ). Applying this result to the case 
s = shows that the eigenvalue fj, 2 ^ = will be realized in a rather large fraction of the total Hilbert space, 
namely a subspace of dimension (2 AT + 2)!/((JV + l)\(N + 2)!) ~ A N+l / (0F7V 3 / 2 ) of the total 4 JV -dimensional 
space. For any N, the mean value of ^ (treating all states as equally probable) is equal to 3iV/16 > 0. The 
standard deviation of /i| is found to be = jL^/3(2 N 2 + N) ps 0.15 (for large N). In the large A^ limit, 
one would recover the result of the classical treatments of Refs [23l [24] . 

Let us also note that, if one were considering the coset- type mass (|4.39[) . including the gravitational- root 
contributions, the result ()4.23|) shows that one would obtain 

~ 2 (ix) 3 i-r 5 3 i + r 5 3 

Mcosot 8 2 8 2 8 ' 1 ' 

i.e. a constant, c-number result, all over the Hilbert space. This result is reminiscent of the full coset result 
(|4.31[) . whose r.h.s., is the quadratic Casimir invariant of the spinor representation of K(G), which is a c- 
number in any irreducible representation. Note, however, that, though the representation space of the spinor 
representation of the Kac-Moody maximally compact subgroup K(G) is finite dimensional, there is an infinite 
number of roots a, and a corresponding infinite number of generators J*, that one should sum over in the coset 
result Eq. (|4. 31 [) . so that the value of /2 2 osct is formally infinite, or, at least, ill defined as an operator. This 
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indicates the need to renormalize it by some ordering prescription. This suggests that, in spite of contrary 
appearances, the renormalized value /i 2 osct = is an allowed possibility. 

Let us come back to the result (|4.43j) predicted by the usual Einstein-Dirac action, and further discuss its 
meaning. Let us first note that the non-zero value (|4.43p of /2 2 is a quantum effect, which directly comes from 
the anti-commutation relation (|4.7I) of the quantized spinor x- If we reinstate the h on the r.h.s. of (|4.7p . we 
see that the objects L Q that satisfy the unit-normalized Clifford algebra (|4.10p are actually related to x v ^ a: 
Xa — 5 VhT a . As p 2 is quartic in \ this shows that /if; is actually | 2 5 h 2 , i.e. of order h 2 . In turn, this 
shows that, when one is in a quasi-classical (WKB-type) regime where the gravitational momenta are large 
compared to h, the mass term will have only a sub-leading effect. This is consistent with the finding above that 
/j, 2 has an intrinsic quantum ambiguity due to ordering problems: indeed, the constant c in Eq. (|4.36j) is also 
seen to be 0(h 2 ). On the other hand, when considering the Klein-Gordon-type equation defined by (|4.36|) with 
Tip = —ihd/d(3 a , we see that a same factor h 2 appears in front of the /3-space d'Alembertian — G ab dp a dp b and 
of both the ^-generated mass /if; , and any quantum-ordering contribution 4c. This means that the mass term 
/t 2 will be important when the wave function <£>(/?, <r) has a characteristic scale of variation in /3-space of order 
unity. 

Let us finally note that the ambiguity in the Einstein-Dirac action linked to using a first-order or second- 
order formalism is equivalent, in view of Eq. (|4.I9[) , to adding to /2 2 a term proportional to —/if; — ^Agrav) 
where /i 2 rav is the contribution to the coset result (|4.3ip coming from the three gravitational roots, i.e. the 
second term on the r.h.s. of Eq. (|4.39l) (which is the sum /if; + /2g rav ). 

This shows that the addition of a suitable multiple of the extra contribution (|4. 19[) to the Einstein-Dirac 
action can modify the basic result /if; in a more general result of the type 

(l-fc)/i| + -^ rav + 4c, (4.45) 

where we took also into account a possible (c-number) ordering constant. For instance, if k = j this would 
generate a mass term | (/2| + /tg rav ) + 4c. When c = 0, this would correspond (modulo the factor |) with the 
coset-like result (I4.44|) . On the other hand, if we choose k = | and 4c = —3/32, we end up with a vanishing 
total squared-mass. 

Keeping in mind all these (quantum) ambiguities in the value of ft 2 , we shall now discuss in more detail the 
quantum dynamics of the Bianchi-IX-spinor system. 



4.6 Quantum dynamics of the Bianchi-IX-spinor system 

The Wheeler-DeWitt equation for the Bianchi-IX-Majorana-spinor system has the form 

\ G ab d 0a d Pb + V^{P) + X> cclltrif (/3, x) + v^ av (p, x) + 

+V smass (f3, x) + \ Mt 2 ot(x)) HP, <r)=0. (4.46) 

Here, the spin- independent term V^ IX (/3) is the usual bosonic type-IX potential (|2.3ip . while there appears three 
different spin-dependent potentials. Using our results above, their explicit expressions are (for type IX) 

i>ix ntrif (/?,x) ' ( ' ' ! 



32 Vsinb^CSi - p 2 ) sinh"(/3 2 -/3 3 ) sinh^(/3 3 - ft) 

i-r 5 

2 ' 

>ix <o ,,\ _ 1 e„-2Si , -2/3 2 i „-2&jv„--p -p ^ 1 + 



(4.47) 



VII^X) - -j( e - 2 ^+e-^ + e -^)( iri r 2 )^, (4.48) 
v£Ls(M = +me-C*+fc+*>(ir 2 r 3 )±±^. (4.49) 
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The various possible values of fi% ot have been discussed in the previous subsection. 

The discrete spin label <r, indicated as argument of takes four values, and reminds us that the various spin 
operators T5, ir^Ta and 1T2T3 act upon a 4-dimensional Hilbert space. In other words, Eq. (|4.46|) is similar to 
the second-order form of the Dirac equation coupled to an external electromagnetic field, namely 



[idu-e A b 



(4.50) 



where er MI/ = | is the spin generator (in the Clifford algebra) . 

In the case of the usual Dirac equation (|4.50|) , the spin coupling term | e er Miy F^ v (linked to the magnetic 
moment of the electron) couples the 4-different components of the (first-quantized) Dirac spinor ip, and embodies 
physical phenomena such as the precession of the electron spin in an external magnetic field. Similarly, we can 
think of the wave function of the universe <fr(/3) as a column vector of 4 components $ CT (/3), which propagate 
in the Lorentzian /3-space, and are deflected by the spin- independent potential Vg X (f3), together with the three 
spin-dependent potentials (I4.47[) - (|4.49[) which are analogous to the \ e a^ v F M „ spin-coupling term in Eq. (|4.50p . 
In addition, the mass term ^tot(x) is a l so spin-dependent, as discussed in the previous subsection. 

In this work, we shall focus on understanding the dynamics of the universe's multi-component wave function 
$(/3) in the BKL regime where one approaches the singularity, i.e. in the limit where the various /3 Q 's all tend 
to +00. For concreteness, let us first consider the case where ji 2 is given by the Bianchi-IX prediction (|4.43|) 
together with a possible ordering constant, but without any further contribution of the type of Eq. (|4.19[) . In 
that case we have 

. 2 3 1 - r 5 



4c. 



(4.51) 



Let us start by noting that the helicity operator T$ commutes with all the terms in the Wheeler-DeWitt 
equation (|4.46l) . This means that the 4 components of $ can be viewed as the superposition of two independent, 
2-component wave functions; say $+ and with 



i±r. 



$ . $ = $ , 



and where and $_ undergo two uncoupled evolutions, 
eigenstate T 5 $ + = +$+) satisfies (with Up = G ab dp a dp b ) 



More precisely 
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(which is a positive-helicity 
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while $_(r 5 $_ = -$_) satisfies 



□/» + K IX (/3) + 
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$_(/?) =0. 



(4.54) 



These two sub-equations have rather different structures: (i) they contain a different mass term, (ii) the 4>_ 
equation is spin-independent (i.e. the two independent components of $_ satisfy the same equation), while the 
two-independent components of $+ are coupled via the presence of the iT\T2 and il^I^ spin-coupling terms; 
and (iii) they contain different spin-averaged potentials (besides the mass term) . Concerning the last property, 



the spin-averaged potential term for 3>_| 
i r a rf,, a 7^ b, are ± 1) is 



(taking into account that (iT a Ti,) = + 1, i.e. that the eigenvalues of 
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while the spin- averaged potential for $_ is 



V-{p) = V*(f3) + - 



32 Vsinh 2 (ft - ft) sinh 2 (ft - ft) 
+ -T^ Tt") • ( 4 - 56 ) 



sinh 2 (ft - ft) 

We note that the spin-averaged potential V+(/3) for $ + is the usual type-IX potential for a diagonal metric, 
which contains only gravitational walls, but no symmetry (or centrifugal) walls. Classically, this potential 
(approximately) confines the motion of the /3 particle (when considering the near-singularity limit) to stay 
within the Lorentzian billiard chamber defined by the condition that the three gravitational wall forms be 
positive: 

tafaaCS) =2 A >0, w a 231 (j3) = 2 ft > , u>g ia (/9) = 2 ft > . (4.57) 

At the quantum level, this confinement mechanism will be blurred by tunnelling effects. As is well-known 
the physics of tunnelling in quantum cosmology crucially depends on the choice of boundary conditions in 
configuration space, see, e.g., [37] ■ Leaving a discussion of more general states to future work, we shall focus 
here on wave functions which are quantum analogs of the classical billiard motions within the appropriate 
billiard chamber (e.g. the chamber (|4.57|) when considering the $+ component). For such states, taking the 
ordering constant c in (|4.53p to its naive value, c = 0, we can qualitatively describe the quantum evolution 
of $ + by looking at the various potential terms in (|4.53p . First, we note that, when trying to penetrate the 
gravitational walls (I4.57|) . i.e. when exploring regions where some of the wall forms, e.g., 2 ft, become negative, 
the spin-averaged bosonic potential V+(/3) = (ft will grow like + | e _4/31 , and will therefore dominate over 
the corresponding (non positive-definite) growing gravitational spin-dependent potential cx e _2 ^ 1 ir 1 r 2 . In 
addition, the spin-dependent potential related to the mass m of the spinor becomes exponentially small in the 
near-singularity limit where the volume of the universe, cx y/g — e~^ 1+ ^ 2+ ^ 3 \ tends to zero (i.e. when the 
sum ft + ft + ft — > + oo). This discussion shows that there will exist quasi-classical spinor- like wave functions 
$+(/?) (with two independent components) consisting of WKB-likc solutions approximately bouncing between 
the walls of the chamber ()4.57|) . and decaying in the "forbidden" domain, ft < 0, ft < 0, ft < 0. Compared to 
the usual, pure gravity Wheeler-DeWitt equation for Bianchi-IX, which would be 

1 

4 



□/j + O 8 ) $oG8) = 0, (4.58) 



for a scalar- valued (one component) <f>o(P), the interesting new feature of the gravity-spinor system is the 
presence of additional spin-dependent couplings. In the case of the equation (|4.53[) . these are the terms 
containing i T1T2 and i r 2 T3. As these two (Clifford algebra) operators do not commute among themselves, their 
presence means that the two independent components of $+ will continuously mix under the influence of these 
terms. We thereby end up with the picture of a quantum fermionic billiard where the various polarization state of 
a spinor-valued wave function mix as the quantum state bounces within the billiard chamber. We shall discuss 
below the precise link between such a quantum fermionic billiard, and its Grassmann-valued correspondent, 
studied in Ref. [22"] . 

In the BKL limit ft + ft + ft — > +00, the term proportional to m in the evolution equations for $ become 
negligible. Let us then consider the special case where m = 0, which should also describe the general behaviour 
when ft+ft +ft — > +00. For the case m — (and c = 0), the evolution equation (|4. 53[) for $ + further simplifies 
in that it contains only one spin-dependent operator, namely iTiIV In that case, the idempotent operator 
iFir 2 commutes with the Hamiltonian, so that we can further decompose into two (scalar) components, 
say 

$V ±) := 1± ^* + , (4-59) 
that evolve independently of each other. Namely, these components satisfy 

(- \ Up + V^(P) T \ (e^ + e-^ + e"^)) $?»(/?) = . (4.60) 
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In other words, when m = 0, we can reduce the dynamics of the multi-component wave function $(/?, a) to the 
uncoupled dynamics of several separate components satisfying different scalar, Wheeler-DeWitt equations. 

Turning to the dynamics of Eq. (I4.54[) , we already noted that it did not contain spin-dependent terms. 
This means that, similarly to the m — , <I> + case discussed above, the two components of $_ satisfy (even 
when m ^ 0) the same scalar, Wheeler-DeWitt equation (I4.54[) . The latter equation has two interesting new 
features compared to the $+ one. First, though the off-diagonal, gauge-like Euler-angle degrees of freedom 
have been eliminated by the diffeomorphism constraints, they have left behind new centrifugal contributions, 
oc 1/ sink 2 (fa — fa) + cyclic, to the potential (see Eq. (|4.56l) . These centrifugal terms create infinite potential 
barriers located at the three symmetry walls: fa = fa, fa = fa and fa = fa. To examine the effect of these 
barriers on the quantum wave function, let us consider, say, what happens near the fa = fa symmetry wall. 
Locally, we can change coordinates in /3-space and use x = fa — fa as new coordinate. As the hyperplane 
fa = fa is timelike (i.e. its normal is spacelike), the coordinate x is a spatial coordinate in the Lorentzian 
/3-space. This means that the $_ equation has, near the x = hyperplane, the structure 



where we completed the spatial coordinate x by another spatial one, y, together with some time-like variable 
t, and where a 2 denotes a positive numerical constant linked to the coefficient of the centrifugal terms. [In 
the case considered here we have a 2 = 1/16, but it is instructive to leave its value arbitrary] Here, t,x,y are 
some linear combinations of the three variables fa, fa, fa, chosen so that the Lorentzian metric G a bd/3 a dfa = 
—dt 2 + dx 2 +dy 2 . We can separate the motion w.r.t. t and y, i.e. $_ = e lhy ~ %ut ip(x). Finally, the behaviour of 
(f(x) near the singular point x = is given by d 2 fix) ~ a 2 x~ 2 f(x). This implies that <p(x) ~ X s with a power 

s satisfying the indicial equation s(s — 1) = a 2 , so that s± = ^ ± \J \ + ot 2 . The important point is that the two 
possible solutions satisfy s + > 1 and s_ < 0. If we reject the possibility of a singular wave function at x = 0, 
this eliminates (x) ~ x s ~ . Finally, we conclude that the centrifugal terms force the wave function to vanish 
at each symmetry wall like \fa — fa\ s+ , \fa — fa\ s+ and \fa — fa\ s+ , with s+ > 1. In other words, the centrifugal 
terms confine the evolution of $_ (fa to six separate chambers, corresponding to the six possible orderings of 
fa, fa, fa: namely, fa < fa < fa, fa < fa < fa, etc. In each symmetry chamber, the WKB-like solutions for 
$_ will bounce between two "hard" symmetry walls (e.g. fa = fa < fa and fa < fa = fa for the chamber 
fa < fa < fa) and the "soft" (exponential) gravitational wall (e.g. the 2 fa — wall for the fa < fa < fa 
chamber) . 

If we use the naive ordering constant c — (which led to a "massless" Wheeler-DeWitt equation for $ + ) 
we see in Eq. (|4.54[) that the $_ components have acquired a positive squared-mass ^i 2 _ = 3/32. This mass 
term will strongly affect the near-singularity behaviour of the $_ wave function. Indeed, as there exist wave 
packets of the free massive Klein-Gordon equation (Dp — = that are approximately localized (in /3 

space) around, say, a time-like line fa — v a r + fa^\ with G a bV a v b < 0, such wave packets yield approximate 
asymptotic solutions of the Wheeler-DeWitt equation (|4.54|) . because, as fa + fa + fa — > +oo (within, say, 
the chamber fa < fa < fa) the potential terms will become negligible in the domain where the wave packet 
is approximately localized. This is the quantum version of the result of [531 121] that a positive /x 2 ultimately 
quenches the chaotic billiard motion, in the near singularity limit . 

However, let us note that, when c = 0, the fact that the q-number /x 2 admits the eigenvalue zero means that 
there will always exist a part of the quantum wavefunction (namely the one described by its $+ component) 
which will exhibit a chaotic behavior near the singularity. In other words, at the quantum level, the classical 
result of [23l [24] does not prevent part of the quantum reality to behave chaotically. 

5 Quantum versus Grassmannian fermionic billiards 

The two crucial new features brought by coupling gravity to a spinor are: (i) the appearance of new spin- 
dependent potential terms; and (ii) the possible presence of a squared-mass term in the Klein-Gordon- Wheeler- 
DeWitt equation. As discussed above, though the feature (ii) is important for knowing whether the mixmaster 




(4.61) 
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chaos is affected or not by the coupling to a quantum spinor, it is also a delicate quantum effect, which is 
sensitive to the choice of ordering in the Wheeler-DeWitt equation, as well as to other choices. We shall assume 
in this Section either that we have made choices so that the relevant squared-mass term vanishes^!, or that 
we are considering situations where one can distinguish the effect of the spin-dependent potential from other 
effects. Under this assumption, we can focus on the effect of the spin-dependent potential terms, and contrast 
a quantum treatment of these spin-dependent terms, to a treatment where the spinor variables are considered 
as Grassmann-valued (G-numbers). 

The recent work [35] has studied in detail the fermionic billiards defined by the chaotic dynamics of G- valued 
(spin-3/2) spinor fields, as they arise in the near-singularity limit of supergravity (both in dimension eleven, and 
dimension four). The two main results of Ref. [32] that are relevant to our present study are: (1) the spin-3/2 
billiard defined by supergravity factorizes into a spin-1 (vector) billiard, and a spin-f/2 (Majorana-spinor) one; 
in view of the results of Ref. [24], it is the latter spin- 1/2 billiard which is of relevance for us here; and (2) the 
spin-1/2 billiard consists of a succession of generalized Weyl reflections, defined in a kind of spin-covering W s of 
the standard Weyl group of the Kac-Moody algebra associated with the considered supergravity model. More 
precisely, it was found that each collision of the universe on a wall labelled by a Kac-Moody root a (i.e. so that 
the corresponding wall form is simply w{0) = a(/3)), corresponding to a spin-dependent potential proportional 
to exp(— a((3)) i x^ Ja Xi causes the G- valued spin-1/2 field x to "rotate", in the spin-1/2 representation space, 
by an angle of ir/4 along the axis defined by J s a . More precisely, the value of the Dirac spinor x. after the 
collision on the wall a differs from its incident value by the matrix 

K a =e £ <*% J «. (5.62) 

Here, e a = ±1 is a sign which will be defined below, and the object J^, which was introduced above, is an 
anti-hermitian generator which represents (in the spinor representation s) the "rotation" E a — E- a of the 
maximally compact subgroup K of G associated to the root a. For instance, in the (3 + l)-dimensional case of 
relevance here, the generators J* associated, respectively, to the symmetry walls and the leading gravitational 
walls, were written down in Eqs. (I4.37[) . (|4.38[) . Note that the generators all include a factor 1/2, so that 
their eigenvalues are ±i. [In conjunction with the ^ prefactor, this implies, as stated, a rotation angle of ±f .] 
Note also that all the generators J* belong to the Clifford algebra defined by the usual (3 + l)-dimensional 
gamma matrices, so that the spinor-wall-reflection matrices (|5.62j) are 4x4 matrices acting on the usual Dirac 
spinor space. [More precisely, they are real, and act on the real Majorana spinor space.] 

The picture of fermionic billiards emerging from Ref. |22j is a growing, chaotic succession of spinor reflections 

T^-Ci n *^^CE«_l ' ' ' ^^Q2 (5.63) 

acting on the spinor index of the G-valued spinor x- It was found in Ref. |22) that the multiplicative group 
W s defined by spin-1/2 billiard, i.e. by all the matrix products (|5.63|) . is of finite cardinality, both in the eleven- 
dimensional case, and the four-dimensional one. For example, in the 4-dimensional case of relevance here, the 
group of products (j5.63[) of the basic reflection generators, Eqs. (|4.37[) . (|4.38[) . associated with the symmetry 
and gravitational walls is a finite group of cardinality 4 x 48. Actually, it happens that the 4-dimensional case 
is quite special in that the symmetry generators happen to commute with the gravitational wall ones (which, 
themselves, reduce to only one element as the three generators (|4.38|) differ only by a cyclic permutation of 
123, which does not change the value of 7 123 ). Due to this, the group W s is the direct product of two separate 
groups. 

Having recalled the results concerning the near-singularity limit of the dynamics of G-valued spinors coupled 
to a chaotic bosonic cosmology, we wish now to clarify the connection of these results to the case of quantum 
spinors, sourcing a cosmological model. In that case, we must replace the classical evolution of a G-valued 
spinor x, considered along the classical evolution of the bosonic billiard described by the /3-particle bouncing 
between symmetry and gravitational walls, by a solution of the coupled multi-component Wheeler-DeWitt 
equation. However, in order to be able to compare the two treatments we must consider a case where the 
Wheeler-DeWitt equation do correspond to the spinor coupling terms studied in Refs. [24] [22]. Indeed, the 
latter works considered classical Hamiltonians of the form 

12 This would, for instance, be the case for the <E>+ components of <E>, with the naive ordering constant c = 0, see Eq. 114.531 1. 
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H = t G ab n 0a 7T ft + J2 ^- 2aW) nl + J2 ^- Q («H Q S Q + c tot (5.64) 

a a 

Here, the n a 's are dynamical variables whose meaning is different if one considers a gravity model or a coset 
model, the E a 's are the spin-coupling terms corresponding to a general root, namely 

V ol = i X ] J s a X (5-65) 

where the J^'s are the products of usual gamma matrices defined in Eqs. (|4.37[) . (|4.38|) for symmetry and 
gravitational roots, and where the constant ctot (defining the squared mass /i 2 = 4ctot) includes various types 
of contributions, which depend on the model considered (as discussed above). 

The important feature of the spin-couplings included in (|5.64[) is that the same quantity n a enters the 
bosonic potential wall ie~ 2Q ^n 2 , and the corresponding spin-dependent potential term ^e~ a ^H a T la . One 
can verify that this is indeed the case for the general Hamiltonian (|3.72l) . both for the terms linked to the 
dominant gravitational walls (with, e.g., 0:123 = 2/3i, n ai23 = — C\ 3 , and </ Ql23 given by (j4~38f ), and to the 
symmetry walls. However, in the latter case, as Eq. (|5.64j) refers to an Iwasawa decomposition of the metric 
one must reconsider the far- wall limit of the Gauss-decomposition-based Hamiltonian (|3.72j) . In that limit (e.g. 
when fa — Pi ^> 1) one finds that the various hyperbolic functions in Eq. (13.79)) do yield a structure compatible 
with the general Iwasawa result (|5.64|) (with, e.g., n Ql2 = ir w n). 

In order to see more clearly the differences between the Grassmannian treatment of spinor couplings, and its 
quantum treatment, let us consider the quantization of the simplest version of the general Hamiltonian (|5.64p . 
namely the case where there is only one wall, corresponding to one root a. In that case, the dynamical variable 
Tl a is clearly seen to be a constant of motion, and can therefore be considered as being a c-number both in a 
classical and a quantum treatment. [In the quantum treatment, we consider eigenstates of the operator Tl a .] 
It then remains to quantize /3 and X- Canonical quantization of the /3 dynamics yields = —id/3, while the 
quantization of \ is done according to the anticommutation relations (|4.5|) , or (|4.7|) in the Majorana case that 
we shall consider here. This yields a quantum Hamiltonian of the form 

H = - J G ab d Pa d Pb + V 2Q («n2 + ie- Q («n Q £ Q + ctot (5.66) 

where S Q is the quantum operator defined by replacing \ in Eq. (I5.65|) by its quantized version, submitted 
to (fC7) . 

As explained above, the quantum wavefunction $(/3, <r) has both continuous indices related to the j3 dy- 
namics, and a discrete one, a, related to the Clifford algebra (|4.7j) satisfied by the quantum \- Let us consider 
solutions of if$(/3, cr) = describing wavepackets colliding on the a potential wall. One can write such solu- 
tions explicitly by separating the variables in Eq. (|5.66|) . Namely, as the potential terms depend only on the 
combination a(f3) of the /3's, we can look for wavefunctions $ of the form 

$(M = e ik ^F(/3 x ,<r) (5.67) 

where the linear form fell describes the /3-space momentum parallel to the wall a(j3) = (i.e. G ah k\otb = 0), 
while F(f3j_,o~), where f3j_ := a(/3), describes the motion perpendicular to the a wall (as well as the spin 
degrees of freedom). Inserting this separated wavefunction in the Wheeler-DeWitt equation H$>({3,<r) = (and 
considering a root a(/3) of norm G ab a a ab = 2, as is the case for all the roots considered here) yields the following 
equation for F: 

d} x F(p u a) = (e~ 2 ^n 2 + e-^n a E a - Q)F(fi ±) a) (5.68) 

where Q is a separation constant which involves both (fe") 2 := G ab feafejJ and the constant c to t in Eq. (|5.66p . 
If S Q were a c-number, this would be the Schrodinger equation in a Morse potential, which is well-known to 
be reducible to the general confluent hypergeometric equation. In our case, S Q is an operator in spin space 
(i.e. acting on the index a labelling the four components of the wavefunction $). However, we can reduce our 
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problem to a c- number- valued S a by first considering wavefunctions $(<x) that are eigenstates of this operator: 
say E q $(<t) = S Q $(d). We can also simplify the resulting equation by shifting j3± by a constant so as to 
absorb the II 2 factor multiplying e~ 2l3± . Namely, one defines /3' ± such that e -2 ^ = II 2 e~ 2 ^ This makes also 
IT a disappear from the second term, except for its sign: e a :— sgn(II Q ) = ±1. When Q is positive, there is a 
solution of the resulting equation 

dj, ± F(p' ± ,a) = (e" 2 ^ + e-^e a E a - Q)F(f3' ± ,a) (5.69) 

which starts as an incident plane wave F(/3' ± , a) oc exp —iy/Qf3', faraway from the wall, i.e. when /3± ^> +1, 
reflects on the wall located around f3' ± = (with an evanescent wave in the 'forbidden' region f3' ± < 0), and 
ends up, with some dephasing, as an outgoing plane wave F((3' ± ,a) oc exp+i^/Q/3' ± when (3'j_ — > +oo. Using 
the notation U(a,b,z) for the second Kummer function (see [38], chapter 13), the exact solution of Eq. (|5 .69[) 
describing this scattering process on the combined bosonic-f spinorial wall a(/3) is given by 

F[(3' ± ] = exp[-e- "i] e" 1 + i £ce I] Q + i VQ, 1 + 2 i VQ, 2 e~^] (5.70) 

We can then extract from this exact expression (using the expansion of the Kummer U(a, b, z) function near 
z = 0) the phase factor e lSa between the incident wave and the outgoing one. We find 

r(| + i £a s Q + i Vq) r(| - 2 i Vq) 

Let us now compare this quantum dephasing with the Grassmannian description, recalled above, of the 
reflection of the spinor Xm on the wall a(/3). [Here, we denote by m, n, . . . spinor indices, to avoid the confusion 
with the use of a as a label for the root] . This description was that the value x' °f t ne spinor after the interaction 
with the wall a is obtained from its value x before the interaction via the following matrix transformation: 

X ' m = (n a ) mnXn = \e £ ^ J «] X n (5.72) 
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where e a := sgn(II Q ) = ±1 as above. 

Similarly to what happens for the usual Dirac equation, we can consider that the various components, 
labelled by to, of the spinor Xm encode the classical polarization state of x- This polarization state can then 
be also encoded by decomposing the spinor Xm onto a basis of eigenstates of the first-quantized hermitian spin 
generator i J*. [We recall that i is i times the product of an even number of ordinary gamma matrices; e.g. 
for the symmetry root an =02 — 01 it is the ordinary first-quantized spin generator i ^7 12 , whose eigenvalues 
are ±i, and whose eigen-spinors are often used to decompose a general spinor into various spin states.] We then 
see that Eq. (15.72)) is saying that the classical dephasing, upon reflection on the wall a, of a classical spinor % m 
polarized so as to be an eigenstate of iJ*, with eigenvalue a a = ±i is given by the phase factor exp— ie a ^a a . 

The quantum dephasing (|5.7ip looks a priori quite different from the classical dephasing exp — ie Q -|cr Q . Let 
us, however, check that they agree, as they should, in the quasi-classical, WKB limit. This limit corresponds 
to considering a large value of a/Q, i.e. high-frequency wavepackets ~ exp ±i\/~Qf3' ± . In addition, we need to 
decompose the quantum phase S a (T, a ) in (|5.71[) into two separate contributions: (i) a spin-independent part, 
S a (0), which can be mathematically defined by replacing S Q by zero on the r.h.s. of Eq. (|5.71[) . and (ii) the 
spin-dependent contribution, 5 a (T, a ) — 6 a (0). This yields a spin-dependent dephasing factor given by 

e i(*«( Sa )- M 0)) = L(l + l £a I] a -zVg) T(±+iVQ) 

r(i-i>/3) r(± + ±£ Q £ a + iv / Q) 

Using now the fact that, for large values of z, one has T(z + a)/T(z) = z a (l + 0(1/ z)), we see that, modulo 
fractional corrections of order l/-\/Q, the spin-dependent part of the quantum phase, i.e. 5 a (T, a ) — <5 Q (0), is equal 
to — ^Sa'Ea, in perfect agreement with the classical result — ^e a a a . This result also shows that the generalized 
spin operators E Q = ix* JaX are the second-quantized versions of the first-quantized gamma-algebra generators 
iJ*, and that classical spinors Xm that are eigenstates of iJ% do correspond, after quantizing x according to 
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Eq. (|4.5p . to quantum eigenstates $ of E a . Note also that the fact that the quantum phase (15.71[) agrees, in 
the quasi-classical limit, with the result obtained by integrating the classical equation of motion of x.i namely 



dtX 



1 



-a(/3)TT 7 



a X 



finally rests on the fact that, according to the general result (|4.8j) . the Heisenberg-picture quantum operator x 
satisfies the same equation of motion. 

We have compared here the quantum and classical descriptions of the 'collision' of the universe on a single 
(bosonic+spinorial) potential wall a{j3). We could discuss, exactly along the same lines, the dynamics of the 
two $+ components of the wavefunction $ in the Bianchi-IX case. Indeed, $+ satisfies the (two-component) 
equation (|4.53|) which (in the case m — 0) is closely similar to the Wheeler-DeWitt equation associated to Eq. 
(|5.66|) . except that we now have spin-dependent collisions on three different gravitational walls. Note, however, 
that the spin evolution is very simple in this case, as the three gravitational- wall spin operators </q 123 , ^q 231 j Ja 312 
all coincide. The spin eigenstates being the same for the three different walls, there is no room for any interesting 
chaotic behaviour of the spin polarization, as could happen if one had had three different (and non commuting) 
spin operators. As for the two remaining components of the Bianchi-IX $ they satisfy the spin-independent 
equation (|4.54j) so that we cannot compare its dynamics to the discussion of [22 . 

There is, however, another Bianchi model for which we can compare the classical and quantum dynamics, it 
is the Bianchi-II model, defined by structure constants C a be — £bcd,n ad with n ab = diag(l,0,0). The classical 
dynamics of this model is discussed in Appendix C, using an Iwasawa fixing of the dreibein. The diagonal 
metric variables are denoted /?i , /?2 , /?3 , while the Iwasawa off-diagonal degrees of freedom are denoted ^12 , ^23 
and V13 . [Note that indices are not naturally cyclically permuted when working in an Iwasawa representation.] 
The quantization of x is done as above, while the quantization of the metric degrees of freedom is conveniently 
done in a j3, v representation, i.e. with momenta conjugated to the three /3's defined as vv a = —idp a , and with 
momenta conjugated to the three i/'s defined as w ah = —id Vab . 

In the Bianchi-II model, there are only two non trivial diffeomorphism constraints (classically given by Eqs. 
(|C.2jl ). in addition to the Hamiltonian constraints (Eq. (|C.I[) ). The quantum mechanical versions of these two 
constraints will then be the following constraints on the quantum state <I>: 



m 12 ® = w 



13 



$ = 



(5.74) 

In the /3, v representation, i.e. for a wavefunction $(/?, v : a) (where, as above, the discrete index a labels the 
spin degrees of freedom), the constraints (|5.74|) ) imply that $(/?, v, a) does not depend on V12 and ^13, but only 
(besides a) on /3±, fii-, and ^23- 

The remaining constraint is the Hamiltonian one, of the form 



H$ = 



(5.75) 



Assuming for simplicity a vanishing Dirac mass, m = 0, but allowing, as above, for an ordering constant c, 
the Hamiltonian operator [whose classical version is (Eq. (jC.Ij) ) ] may be written as 



n 



1 



-4/3i 



+ - [S^l 2 + £M 2 + (E^l + 2 e^-^ro 23 ) 2 
8 L 

+c 



-2/3i 



«X f 7o7i23X 



(5.76) 



It is amusing to note that the quantum dynamics of this Bianchi-II model is simpler to analyze than its 
classical dynamics. Indeed, as discussed in Appendix C, it is difficult to solve generically the classical (or 
Grassmannian) spinor evolution equation (IC.44I) . By contrast, in the Schrddinger- Wheeler-DeWitt picture of 
Eq. (|5.75p there are no fermionic equations to solve. In addition, we have that: (I) the spin dependence of the 
Hamiltonian is analog to that of a symmetric top, i.e. involving the spin operators E ab only through E^ 23 ], and 
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the combination E^ + E^ ; and (2) the additional spin-dependent term linked to the gravitational wall, 
i.e. «X^7o7i23X, commutes with the E ah -dependent terms (see the subsection l4.4[) . [Note that this commutation 
property (2) is valid in the case where one would consider a complex, Dirac spinor \, because of the independence 
of the two real Majorana components of such a general %.] 

In other words, we can diagonalize the Hamiltonian by imposing that the state $ is a simultaneous eigenstate 
of the following three operators: 

E"l 2Sl $ = CT 2 3$, 
/g[i2]2 + g [ i 3 ]2 + g [ 23]2\ $ = 52$5 



^X f 7o7i23X$ = (5-77) 



On each a component of such a state, the Hamiltonian constraint equation leads to an equation of the form 

$ = 0. (5.78) 



By using the results of subsection 14.41 the allowed values of the quantum numbers S23, C g and S 2 , and their 
multiplicity, for a Majorana spinor are, 

[t 23 , C g , 5 2 ] mult . : [0, ±1/4,0]!, [±1/2,0,3/4]!. 

while, for a Dirac spinor they are : 

[(723, C g , S 2 } mult . : [0,0,01a, [0, ±1/2,0], [±1/2, ±1/4, 3/4] 2 , [0,0,2], [±1,0,2]. 

Note that in the Majorana case 023 and C g cannot be both non zero, and we have the simple link S 2 = 3er 2 3 . 
These links are relaxed in the Dirac-spinor case. 

Eq. (|5.78|) can be solved by separation of variables. Indeed, the wave operator Dp depends on three (3 
variables, while the potential walls entering the equation involve only two combinations of the three /3's, namely 
0^123 (/3) = 2/?i and a 23 (/3) = ^3 — Therefore there exists a linear combination of the three /3's which will be 
'orthogonal' (in the Lorentz-/3-space sense) to the two combinations 0123 (/3) and 023 (/3). It is easy to see that 
ao(/3) := 2/3i + (3 2 +/?3 is such a combination. Actually, one can easily check that the three variables ao, 023, 01123 
define an orthogonal coordinate system in Lorentzian /3-space, with oo being a time- like coordinate, and 023, and 
O123 two space-like coordinates. Finally, a generic solution can always be expressed by superposing separated 
solutions of the form: 

*G8,Hh) = e l P^e lk W 1+ ^+^F 1 (2f3 1 )F 2 (J3 3 - f3 2 ) (5.79) 

where the two functions F\ and F 2 satisfy exactly the same Morse-potential-type Schrodinger equation that we 
encountered above (in the single-wall case) namely: 

Fl'pfr] = ( e - 4ft -2C ff e" 2 ' 3l -Qi)Fi[2/3i] (5.80) 
^-A] = (fe 2 ^-^ +a 23P e^-^ - Q 2 ) F 2 [p 3 - fa] (5.81) 

Here the two separation constants Q\ and Q 2 must satisfy the following mass-shell condition 

~k 2 + Q 1 + Q 2 = -hs 2 -2c=-^ 2 (5.82) 

where /i 2 is the squared-mass of the Wheeler-DeWitt equation. If we look, as above, for wave functions that 
vanish behind the 023, and 0123 walls, we must choose Kummer's f7-type solution for F\ and F 2 . 



34 



When Q\ and Q 2 are positive, the physically relevant solutions are (when assuming p > for definiteness) : 



-Fi [2/3i] = exp[-e- 2 ft ] e" 2 ^^U^-Cg + i l + 2i y/Q u 2e" 2 ft ] 



F2W3 - fo] = exvi-pe-^-^} e -<v^(/»s-ftO 

x C/[i(l + a 23 ) + zv / &, l + 2zV^, 2pe 



-(/3 3 -/32)l 



(5.83) 



(5.84) 



Far from the two walls, these modes propagate as plane waves in all the variables, with /3-space momenta w a 

of the form 

zu a [e G ,e s } = {2(k + e G y/&),k + e s \fQ~2, k-e s 



satisfying the mass-shell condition 



Lt H7 a 137j, = —fl = — — O 



4c 



(5.85) 



The quantities eq and £5 in the momenta w a denote some ±1 signs, that flip upon collisions on the walls. 

As above, we can also compute the phase shifts of these modes as they reflect on a wall. More precisely we 
find that, for given quantum numbers C g , Qi, 023 and Q 2 , the phase shift a g of F± as it reflects on a gravitational 
wall, and the phase shift a s of F2 as it reflects on a symmetry wall, are respectively given by: 



r(i/2-e g 



/ Q7)r(l/2 + 2i v /& 



T(l/2 -C g + i VQi) r (V2 - 2 i VQl) 
r(l/2 + (J23/2 - j VQ^) T(l/2 + 2 i yg^ 
" T(l/2 + a 23 /2 + z V&) T(l/2 - 2 i ^&) 



(5.86) 

(5.87) 
(5.88) 



As above, we can deduce from these results the intrinsic phase shifts due to the spin dependence of the walls 
by subtracting the phase shift of the spin zero mode. 

Let us finally note that the existence of a non-zero mass term /i 2 = £> 2 /2 + 4c in the Wheeler-DeWitt 
equation can lead to an interesting phenomenon (whose classical analogue is discussed in Appendix C) . Indeed, 
a strictly positive /i 2 (e.g. corresponding to c = and S 2 ^ 0) forces the classical trajectory of the wavepacket 
to stay time-like in /3-space, i.e. prevents it to reach the /3-space light-cone. Therefore, such a mass term 
constitutes a kind of potential wall that prevents the wavepacket to reach the light-cone. We can therefore 
think of the dynamics described by the Bianchi-II Hamiltonian above as that of a quantum particle moving in 
a three-dimensional Lorentzian space, and confined by three different walls: the two spacelike walls 0123 and 
a 2 3 (that prevent the particle from going on the negative sides of those spacelike walls) , and a third effective /_i 2 
wall that prevents the particle, after it has bounced on the spacelike walls and aims towards the light-cone, to 
reach the light-cone. These three walls thereby define a kind of waveguide that oblige the particle to move in a 
time-like direction, which is somewhere midway between the spacelike walls and the light-cone. The interesting 
consequence of this waveguide phenomenon is that it can trap the particle in a bound state, confined between 
all these walls. This happens when both Q\ and Q 2 are negative (in view of the mass-shell condition (|5.82|) 
this happens only when /x 2 > 2fc 2 > 0). In that case, one should no longer consider scattering states of the 
Morse-potential equations above, but rather bound states in the Morse potentials. It is well known (see, e.g., 
|39j ) that these bound states occur for the following quantization conditions 



V-Qi = ~n\ - -+C g , 



and 



\/ -0,2 = -Tl2 - X (1 + CT 23) • 



(5.89) 
(5.90) 
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where n\ and ri2 must be natural integers (starting with 0). For instance, in the minimal case where the 
ordering constant is simply the naive value c = 0, there will exist only one such bound state, namely the ground 
state m = n-2 = 0, with C g = —023/2 = 1/4, Qi = Q2 — —1/16, and k — ±1/4. This solution furnishes an 
interesting example of a quantum cosmological ground state associated to a Bianchi-II billiard. Note that we 
have discussed here a state which is bound simultaneously within the two separate Morse-potential equations 
associated to the gravitational and symmetry walls. There can also exist semi-bound states, i.e. states which 
are bound w.r.t, say, the gravitational-wall Morse-potential, but which represent scattering states w.r.t. the 
symmetry-wall potential. 

6 Conclusions 

We have studied the minisuperspace quantization of spatially homogeneous (Bianchi) cosmological universes 
sourced by a Dirac (or Majorana) spinor field. In the main text we used a formulation of the spinor dynamics 
in which the local 50(3, 1) local Lorentz symmetry of the vielbein is gauge-fixed from the start. [Appendix A 
compares this approach to the one where one does not initially fix the vielbein.] In the Bianchi types IX and 
VIII (corresponding to simple homogeneity groups G) we fixed the 50(3) freedom in the dreibein by using the 
existence of a three-dimensional automorphism group of the Lie algebra of G. In the Bianchi type II case, we 
fixed the 50(3) freedom of the dreibein by using an Iwasawa decomposition (which happens to be compatible 
with the automorphism group of the corresponding G). 

In the Bianchi types IX and VIII, the quantum version of the three diffeomorphism constraints means that 
the wavefunction does not depend on the three Euler angles parametrizing the (pseudo-)orthogonal matrix 
5 G SO (3) (for type IX) or S € 50(1,2) (for type VIII), entering the Gauss decomposition g — 5 T diag5 
of the metric g. This is the quantum version of the classical possibility of restricting g to a diagonal form, 
diag = diag(e~ 2ft , e" 2fe , e" 2& ) . 

The quantization of the homogeneous spinor (denoted x after a rescaling) leads to a finite-dimensional 
fermionic Hilbert space, which means that the wavefunction of the universe, which, in the bosonic case, has 
only one component, becomes multi-component in presence of a spinor field. In addition, in the Majorana case, 
the four components of the wavefunction can be identified with the four components of a spinor of an Euclidean 
0(4). 

The multi-component Wheeler-DeWitt equation satisfied by the wavefunction $ is similar to the second-order 
form of the Dirac equation coupled to an external electromagnetic field ( Ui — e A^) 2 + \ e + m 2 ] ifj — 

0), namely it has a structure of the form (for types IX and VIII) 

(- \ G ab dp* d h + V {p) +J2 V « O 8 )^ (X) + E V ™ 8 ) W) 2 + \ mL Ot)) *(M = . (6.91) 

^ a a' ' 

where £ Q (x) = ix JaXi an d ^a'(x) = *X ^a'X are some bilinears in x> and where a and a' run over some 
sets of linear forms in (3 (or 'roots'). There exist a limit in /? space (the far- wall, or BKL, limit) where all the 
potentials Va(/3) tend exponentially towards zero. [The existence of this limit defines the separation between the 
V}i(/3)-terms and the squared-mass term \ A* 2 t(x)-] The main features of this multi-component Wheeler-DeWitt 
equation are the following. 

The squared-mass term /i 2 ot (x) is a quantum effect, which is of order ti 2 , and which is affected by several 
sorts of quantum ordering ambiguities. We discussed the fact that it was different in the original minisuperspace 
Einstein-Dirac theory compared to the spin- 1/2 Kac-Moody coset proposed in [53]. This suggests that we should 
choose an ordering (and additional terms, such as (|4.19[1 ) such that the total squared-mass term vanishes. On 
the other hand, if we do not do so, but use instead the naive ordering that looks natural in the quasi-Gaussian 
spacetime gauge N = ^Jg, one gets a specific prediction for /i 2 ot (x). One then finds that this quantum (spin- 
dependent) operator admits the eigenvalue zero in a fraction of the total Hilbert space. This ensures that a part 
of the total wavefunction, i.e. a part of the total quantum reality, will formally continue to behave chaotically 
near the singularity, in contrast with the case where the spinor source is treated classically, where /x 2 is a strictly 
positive c-number [23] (see (|4.32p ). 
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We discussed in some detail the physical effects linked to the other terms in the multi-component Wheeler- 
DeWitt equation ()6.91|) . We studied in particular the spin-dependent terms of the form J2 a V a {0)E a (x). Such 
terms appear both in the Bianchi IX and VIII cases, and in the Bianchi II one. In the Bianchi IX and VIII cases 
the set of roots a entering these terms are the three gravitational roots, while, in the Bianchi II case there appear 
both a gravitational root and a symmetry root (see (|5.78|l ). When combining these spin-dependent terms with 
the corresponding potential terms ~ e - 2a (P) of the spinless potential Vo(/3) involving the same root a, we found 
that they lead to a Schrodinger equation in a Morse potential. By studying the quantum scattering on such 
spin-dependent Morse potentials we could relate the quantum dynamics of wavepackets reflecting on them to a 
previous study of fermionic billiards, which used Grassmann-valued spinor fields. The quantum spin dynamics 
of Bianchi IX and Bianchi VIII happens to be rather trivial because all the corresponding spin-dependent 
couplings S Q (x) can be simultaneously diagonalized. A more interesting spin dynamics would, however, be 
obtained in more complicated models (e.g. in higher- dimensions) where the various E a (x)'s do not commute 
among themselves. 

We also studied in detail the Bianchi II model. For this case we could provide the exact general solution 
of the quantum dynamics. It can indeed be decomposed in separated modes, which can all be expressed in 
terms of confluent hypergeometric functions. Some of these solutions describe wavepackets reflecting on the 
gravitational and symmetry wall forms CH123 (/?) and 0:23 (/?) of the Bianchi II model, while other solutions (present 
when /i 2 > 0) can describe interesting bound states, trapped between the walls 0:123 (/?) and 0:23 (/?), and the 
effective wall generated by the positive squared-mass term. 

Note finally that the Appendices provide more details about several formal aspects of the gravity-spinor 
Hamiltonian dynamics, as well as a study of the classical limit of this dynamics. 
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Appendices 

These Appendices are devoted to some aspects of the Hamiltonian formalism applied to the dynamical problem 
constituted by the Einstcin-Dirac equations considered in the framework of homogeneous Bianchi (class A) 
cosmological models. As described in the main text, we use an adapted tetrad constituted, at each point, by 
a time-like vector orthogonal to the slices of homogeneity and a dreibein tangent to these slices: thus these 
dreibeins are defined up to local 50(3) rotations. In the first Appendix we do not fix the local 5*0(3) freedom 
in the dreibeins, and provide all the constraints and Hamiltonian evolution equations of the corresponding 
dynamical variables. In the second Appendix we make use of an Iwasawa decomposition to fix the 50(3) gauge 
freedom and display the Hamiltonian and momentum constraints and general constants of motion, bilinear in 
the spinorial variables. The third Appendix consists of a sketch of the resolution of the classical equations for 
the particular case of the Bianchi II cosmological model. The main aim of this discussion is to provide the 
elements needed to compare the classical and quantum dynamics of a billiard collision near the cosmological 
singularity. 

A Locally S'0(3)-invariant approach to the Einstein-Dirac dynamics 

In the approach |T5l [16] where one fixes the time-like vector of the vielbein but leaves an 50(3) rotational gauge 
symmetry in the choice of the spatial dreibein, the dynamical variables (in an Hamiltonian formalism) are N, 

N a , ht , together with the conjugate momenta to the dreibein and the spinor, say, respectively, H ka , 
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11^4 and U A . Starting from the Lagrangian action L = Leh + Lr>, where 

Leh :=N^g( 3 R + K ab K ab - {K a a f) , 
with K a b denoting the second fundamental form, Eq. (|2.25|) . and 

one computes the two sets of momenta. First the gravitational momenta conjugated to the dreibein 



where we defined 



as well as 



U ka := dL/dh ka = 2 W a +2P ka - S ka (A.l) 



Il ab := +V5 [K ah - g ab K] =Tl ba , (A.2) 

P ab := -V? Wt"^* 4 a/l ? = J V? *7 6 * .9 ab = ^ , (A.3) 

S ab := 1 Vs^V^ = -<S ba , (A.4) 

D* 6 := fc* IP 6 , P* 6 := h\P ab , 5 ib := S rfea /i| . (A.5) 

G G 



Note that, for convenience, we follow Ref. |16j in working in the Appendices with the spinor bilinear S ab , 
which differs by a factor one- half from the one used in the main text: namely S ab = |S ah , with the definition 
(12381) . 

We also have the fermionic momenta that lead to primary constraints : 

U A := d L L/d* A = -^ (* 7 6 ) A = -H^^, (A.6) 

IP 4 := d L L/d¥ A = 0. (A.7) 
The Cartan one-form w = pdq is given by: 

zu = U ka dh ka + T1 A <M A + Tl A d¥ A (A.8) 

and the total Hamiltonian is given by 

H = NH + N k H k + n [ab] H [ab] (A.9) 
where (using when necessary the definition H rs := ft,^n fes ): 

n = (Vn rs - i(rp 2 ) - - + + ^ (A.io) 

and 

H s = -D r n sr + v /5*7 6 £> fl * , -H [afcl = i(^n lcb -/i|n fea ) + s'[ ab ] . (A.ii) 

2 

The coefficients Qr a &] in Eq. IA.9I are Lagrange multipliers measuring the arbitrariness in the local rotation 
rate of the dreibein. We will fix them in Appendix [Bl To the three Lagrange multipliers N, and f2[ ab ] are 
associated the three constraints: 

M»0, H s &0, H [ab] ^0. (A.12) 
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It is useful to use slightly different canonical spinorial variables: \ A := g 1 / i ^! A and cj)\ = ig 1 / i '^\. 
Accordingly we obtain: 

n A d^ A = ^ A d X A - ^ X ~ g dg = 4>W A - \<$X h ka dh ka . (A.I3) 

This redefinition of the spinorial variable entails a corresponding redefinition of the canonical momentum conju- 
gated to the dreibein ht . Instead of H ka the new momentum conjugated to ht becomes 7r fca = LI fca — \<$x h ka ■ 
This shift in the definition of the gravitational momentum has the effect of suppressing the second , 2 P ka , con- 
tribution in the above expression for H ka , so that one finally obtains 

TT ka = 2 U ka -S ka . (A.14) 

G 

When appropriately moving the indices by means of the dreibein /i| or its inverse, the two contributions to the 
redefined gravitational momentum read: 

n ab = + h V a ) > s&b = i^V S x (A.15) 

Let us now apply the above general formalism to homogeneous spaces of Bianchi class A, i.e. with structure 
constants, Eq. (|2.3p . of the form 

C a cd := n ab e bcd . (A.16) 



3 " * ' ab„ ^„2 



In such a framework we have 

^R = -l{n ab n ab ~^n 2 ) (A.17) 

where the indices on n a b have been lowered with g a t,, and where n := g a b n ab . The spatial Dirac derivative of 
has the general form 

^* = ^* + i^ 7 5 75*, (A.18) 

where we have taken into account the possibility to introduce a spatial variation of ^ via constants A a such that 
the Lie derivative of the spinor field along the generators £ a of the isometry group of the homogeneous slices 
verify £^ a ^ = i A a Integrability of these equations implies that A a n ab = (see Ref. [16]). When A Q ^ the 
fermionic degrees of freedom are expressed as a product of a time dependent dynamical factor multiplied by a 
space dependent phase factor. However, this spatial dependence disappears in the stress-energy tensor T^ v (ty). 
As indicated by the explicit i entering the condition £^ a ^! — i X a ^ such generalized homogeneous solutions are 
only possible in the complex, Dirac-spinor case. Seen from the viewpoint where one decomposes a Dirac field 
^> into two real Majorana fields ^r^i, these solutions correspond to having spatial oscillations turning one 
Majorana mode into the other one (but keeping quadratic expressions ~ ^ 2 R + constant in space). 

From all these results we obtain the various terms defining the total Hamiltonian (|A.9|) in the framework of 
Bianchi cosmologies: 

H = -L ( F n u - \ {% f) + ± (n« b n ab - \ri) + i JL fax ~m^ X (A.19) 

where 75 = 70717273 is the common Clifford algebra generator linked to spinorial gravitational walls, see Eqs. 
(|P51l . while 

Us = 2 n™" Crasn + jUV^X GwU + I \ s $ X 
G 4 

= ir^C^sA+iXstfx, (A.20) 
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Here the hatted indices on 7r m ™ and Cmns are obtained from their original forms ir ka and C a bc by suitable uses 
of the dreibein ht or its inverse. 

ka 

The Einstein-Dirac Lagrangian leads to a set of second class constraints that oblige us to consider Dirac 
brackets instead of Poisson brackets. Denoting Dirac brackets as { , }* we obtain from the constraint equations 
(|A.61 IA.7j) the canonical brackets : 

{h fna ^ b Y = 5%5 b a [ x *J B y = 5* = i{ X A ,X^Y (A-22) 
with corresponding classical equations of motion 

ka^iha^Y , TT^pMf}* , ± A « { X A , H}* (A.23) 

completed by the first class constraints: 

H^O « M «0. (A.24) 

So that, in the gauge N a = 0, we obtain (with h := ^J~g and ■k^ := h,%ah^ b ir mb ): 

ka = ^ ((Ha + Hi) (^ b M Ha) + "[fofc] C (A-25) 

from which we deduce 



and 

7T 



To be complete let us also display the Dirac equation : 

X = N (i 7 °A- ^75 - m 7 6 ) x + l%^ [&h] X ■ (A.29) 

B Approach using an Iwasawa representation to fix the local 5*0(3) 
gauge 

In Appendix [3 we displayed an Hamiltonian formalism in which the local rotation rate of the spatial dreibein 
is left unfixed. Starting from this formalism, we can then gauge fix the local rotation of the dreibein by choosing 
a specific way of defining a dreibein ht from the metric components g a b- In the main text, we did this by using 
the form (|2.5I) . with a matrix S belonging to the automorphism group of the Bianchi structure constants C a be- 



— (hn„h q ~ — ht ht) = f2r.fi 

2 \ a 1 b bq aj [ab] 



(ab) jy 



(jlac 



hi + ht hi 

b bq a 



2h c 

— o-t 
N ab 



(A.26) 
(A.27) 



U f^qaib i „qabi ^qbia 

~ Hob — T I 7T 7T + TT 7T — TT H TT 

Ah V 

2m" 6 ) 



(A.28) 
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Here, we show the form of the Hamiltonian obtained when fixing the definition of the dreibein by means of an 
Iwasawa representation, namely 



/ e~* \ / 1 ua vi 3 
ha = e-A- 1 v 23 

\ e-fc ) \ 1 

i.e. 

h-ka = e ~^ kl, ka without summation on fc. (B.2) 

Let us emphasize that there are (only) three (non- vanishing) off-diagonal v variables, namely v\ 2 , ^23, ^13, i-e. 
the Vka' s such that fc < a. [The ^ a 's with fc > a are equal either to or to 1.] These three metric off-diagonal 
variables are the Iwasawa-decomposition analogs of the three Gauss-decomposition Euler-angle variables. The 
momenta conjugated to the Iwasawa diagonal and off-diagonal variables, /3j. and v^ a (fc < a) will be denoted by 
w k and vD ka (fc < a): (zd 1 ,m 2 ,vo 3 ,-cu 12 ,m 13 , to 23 ). Note that w ka is only defined for fc < a. This contrasts with 
the momenta ir ka conjugated to the dreibein which are defined for all indices: fc < a, fc = a and k > a. 
Among these, the ones with fc < a and k — a are easily related to the w ka (fc < o)'s and the f3's by considering 
the Cartan one-form (IA.8I) . Specifically, one finds : 

whenfc <a : Tr ka = e h w ka , 

Hu ( k 1 \ ^ ka \ 

TT = 6 yVJ W V ka) ■ 

k<a 

On the other hand, the calculation of the remaining components of the gravitational momenta, i.e n ka for fc > a 
is less immediate and must be based on the fact that, when viewed within the non- gauge-fixed formalism of 
Appendix A the Iwasawa gauge (|B . 1[) introduces three new constraints: ht — for fc > a, that do not commute 

with the constraints The presence of these new constraints oblige us either to apply the Dirac procedure 

or to solve them explicitly. Using this last option allows one to compute, recursively, the remaining momenta 
ir ka with fc > a, i.e. tt 21 , tt 31 and 7r 32 , by using the relation 

TT ka = h^ d h k d - 2S [ka \ (B.3) 

starting from (fc, fc — 1) down to (2, 1). Indeed, noticing that hp = if p < a and h k d — if fc > d, we see that 
for given values of fc > a, the right hand side of Eq. (|B.3|) involves only already known momenta 7r^ d with either 
p > k or k >p > a but d > p. In three dimensions, this yields the explicit results 

tt 32 = e^-^vu 23 + e^v 23 vo 3 + 2e h S^ 

tt 31 = e~^\e 2 ^(y lz - v l2 v 2Z )vo 3 + e 2 ^w 13 - e 2 ^v 12 w 23 ] + 2[e^S^ - e^v 12 S^] 
tt 21 = e^[*W + e" 2 ^-* 31 V 2 + "2^ 13 ) + (^23^12 - ^i 3 )^ 23 ] + 2e^S^ 

(B.4) 

Now we have in hands all the ingredients needed (|A.14[) to write the remaining Hamiltonian equations and 
constraints (|A.19IA.20p in terms of the Iwasawa variables. 

The Iwasawa Hamiltonian density (of weight +1) reads, when using a diagonal forrrF 3 ! n ab = diag {ni, n 2 , 77.3} 

13 As is well known, the matrix n ab can always be diagonalized; we assume it has been. 
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for the symmetric matrix n entering the structure constants C a d = n ab Ebcd, 

H = e i: i A|i(n7 1 ) 2 + (^ 2 ) 2 + ( W 3 ) 2 -^(^ 1 +^ 2 + ^ 3 ) 2 

(e^-^(zu 12 + v 23 w 13 ) + S^) 2 + (e^-^m 13 + S^) 2 + (e^'V 3 + S^f 



-n\n 2 [e 
-nm 3 [e 



-n 2 n 3 [e 



-2(/3 2 +/3 3 ) 



-4/3, 2 



12J 



-4/3i„2 



13 



-2(/3 1 + j 8 2 ) 2 



23. 



+ e 



-4/3 2 „? 



23 



-4/3i 2 2 
^12^13 



+e -2(ft+fe) (l/ 2 3 _ 4iyi2 ^ ^ + ^ ^ 



-i[n ie - 2 ^ + n 2 (e- 2 * + e^ 1 v 2 2 ) + ^(e" 2 * + e" 2 ^ 2 , + e~ 2 ^v 2 z )]tf l5X } 

while the momentum constraints reduce to 

H 1 = ^ {n 3 [e 2 ^-^m 23 + v 13 vo 12 - v 23 (vo 2 - vo 3 + v 23 vo 23 ) + 2e^~^S^ 
-n 2 (w 23 + v 12 zu 13 )} 



(B.5) 



(B.6) 



n 2 



e 2(ft-/3 3 ) ^13 + e 2(ft-&) ^ (n7 12 + ^ w l 3) 



e^ 2 {ri! n7 13 + n 2 i/i 2 O 23 + ^i 2 n7 13 ) 

— Tl3 

+i/ 13 (t^ 1 - ^ 2 + 1^3 tn 13 + ^ 23 tn 23 ) + 2 e^-^V 23 ^ 12 J - 2 e^"« S^] } (B.7) 
e" 8 {-m [n7 12 + z/ 23 n7 13 ] 

+n 2 [e 2 ( ft -^)(n7 12 + ^ 23 w 13 ) - wlw 1 -w 2 + u 12 {vj 12 + v 23 vj 12 )] + 2e^S^} 
-n 3 [v 13 (e 2 ^-^vo 23 + (tu 1 - vo 2 ) + 2e<*-*> S™) + v\ z {w X2 + v 23 vj 13 ) 

-^23 (e 2 ^-^vj 13 + 2 e^-P* S™)]} (B.8) 

The constraint algebra reduces to 

{H,H a } = {Ha,H b } = -C c ab H c . (B.9) 

Obviously bosonic and fermionic equations decouple in the sense that the bosonic sector only involve 
quadratic expressions in the fermionic variables. In particular, using the previous equations we obtain (un- 
der the simplifying condition A a = 0) : 

(<A t 75X)'=2m0t 7 9 75X (B.10) 
which constitute a hrst integral when the Dirac mass m vanishes (chirality conservation), and 

£12 = 1 (>-/3s ^23 _ e ^S 23 W 13 ) (B.11) 

plus two similar equations obtained by circular permutation of the indices 1, 2 and 3. It is worthwhile to notice 
that these spin equations are independent from the Bianchi type considered. They lead to the conservation law 
of the norm of the spin tensor : 



(5 i2 ) 2 + (S i3 ) 2 + (S 23 ) 2 = § 2 = C st 



(B.12) 
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C Classical dynamics of the Bianchi-II-Dirac system 

To analyze the behaviour of the spinor variables under a collision, it suffices to restrict ourselves to the (almost) 
simplest case: the Bianchi type II homogeneous model. The relevant structure constants are obtained by using 
n 1 = 1, n 2 = n 3 = in Eq. (|A.16I) . An important remark is that (see Ref.|16j. in which various special solutions 
of the Einstein-Dirac system are discussed) the subgroup SInv of unimodular SX(3,R) transformations that 
preserves these structure constants is a five dimensional group whose generators will provide, a priori, five 
constants of motion [see eqs (|C.4HC.6|) ]. 
The Hamiltonian constraint is: 



\ L i i 



1 

r 



-4/3i 



(e^-M(zv 12 + v 23 vj 13 ) + S^) 2 + (e^~^zu 13 + S^) 2 + (e^~^w 23 + S^f 



'75X: 



(C.l) 



while the three diffcomorphism constraints reduce (upon using n 1 = 1, n 2 = n 3 = 0) to only two non-empty 
constraints, namely 

(C.2) 



w ~ U 



,13 







Let us notice that the Hamiltonian does not depend on the variables v\2 and ^13, so that there will be no 
ordering ambiguities when quantizing this system. 

When replacing the constraints (|C.2j) within the Hamiltonian, and using as lapse and shift N — e~^ 1+ ^ 2+ ^ , 
N k — , one ends up with a Hamiltonian of the form 



H = - 

4 



2 6 



i i 
+ 1 [( S [i2])2 + (S [i§]j2 + ( e (0 2 -p 3 ) w 23 + £[231)2 

The non vanishing extra constants, taking into account the momentum constraints, are 
w 23 « p , 

uj 3 — w 2 — 2 1^23 vj 23 ~ C\ , 

v 23 (m 3 ~ vj 2 ) + (e 2 ^~^ - v 2 2Z ) vo 23 + 2 e^-^S^ « C 2 . 
Moreover we also obtain that 

is a constant of motion, as well as the norm of the spinor field 

and, if we assume that the Dirac mass m is zero, 

(</4) X f 75X « C s 

constitutes one more constant of motion. 

As usual in this framework we adopt as lapse and shift functions: 



(C.3) 

(C.4) 
(C.5) 
(C.6) 

(C.7) 
(C.8) 

(C.9) 



N = e -(/3i+/32+As) 



N« = . 



(CIO) 
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The equations of motion reduce to F — N [F, %]. When combined with the energy and momentum constraints, 
and the previous constants of motion, they lead to: 

= +lpe^-^S^ (C.ll) 

£[i3] = _}_ pe (P2-p s ) S m (CU2) 

from which we deduce another constant of motion: (S^ 12 !) 2 + (S^ 13 ') 2 = L 2 , in accordance with eqs (|B.12[) and 
(|C.7p .[ Note the links s 23 = ct 23 /2 and L 2 + s 23 = § 2 = S 2 /4 between these constants of the motion, and the 
other constants introduced above; in the following discussion of the integration of the classical equations of 
motion we denote some constants of integration by the same symbols as their corresponding quantum numbers 
introduced at the end of section [5]] 
We also obtain : 

fa-k = \(w 2 -w 3 ), (C.13) 
zu 3 ~w 2 = 2e iP2 ~ M p{e if)2 - p3) p +s 23 ). (C.14) 

Obviously this system may be integrated by two quadratures. However it will furnish only one new constant of 
integration as from eqs (|C.4IIC.6|) we obtain: 

\{vj 2 - vu 3 ) 2 = P C 2 + \C\ + s 23 2 - (pet*"*) + s 23 ) 2 . (C.15) 

i.e. equivalently 

02 - As) 2 = K 2 - (p eW>-M + s 23 ) 2 (C.16) 

where we have put K. 2 := pC 2 + C\ /4 + s 2 3 2 = Q2 + s 2 3- Let us notice that this equation implies that fi 2 — fiz 
is bounded from above if there is effectively a symmetry wall, i.e. unless p = 0. More precisely the domain of 
variation of (3 2 — (3 3 are the following, depending on the relative values of K,(> 0) and s 23 i.e. of Q 2 : 

Q 2 > [K? > s 2Z 2 ) (C.17) 

ifp S23 >0 : 2 -fte]-oo > H(/C-|«23|)/bl]] 

ifp S23 <0 : / 8 2 - / 8 3 e]-cx),ln[(/C + |s 23 |)/|p|]] 
Q 2 < (K 2 < s 23 2 ) (C.18) 

ifp S23 < : /3 2 -/3 3 e]ln[(| S2 3|-/C)/b|],ln[(| S23 |+/C)/b|]] 

The situation correspondingto integration constants obeying the inequality (|C.17[) corresponds to a situation 
where j3 3 is (almost) alwayQ greater than /3 2 . When K. 2 < s 23 2 , the difference between /3 2 and /3 3 could be of 
constant (positive or negative) sign or oscillate around zero. Once integrated, eqs (|C.11|) and (|C.12[) furnish the 
spin components and S [i3] after a single quadrature while the metric coefficient ^ 23 is directly obtained 
thanks to equations (|C.5[) and (|C 13[) : 

1*3 = -(183 -ft)-! 1 . (C19) 
p 2p 

Another constant of motion is obtained from the equation 

zd 2 +m 3 = -2NH = 0. (C.20) 

We denote it by w 2 + zu 3 = 2k . 



'Always if the ratio (K. =p |s23|)/|p| < 1, otherwise there is a short period of time during which P2 > fls- 
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The remaining bosonic equations, taking into account the previous ones, reduce to : 



02 + 03 

V\2 

*>13 



—w 1 — —k 
4 2 



~2 W ' 



2e" 4ft -2C g e~ 2 ^, 

5 [i3] e (ft-fe) +e (ft-ft) 5 [i2]^ 3 



Equations (jC.21[) and (|C.22|) immediately provide: 

20i + 02 + 3 = k (t - t) . 



(C.21) 

(C.22) 

(C.23) 
(C.24) 
(C.25) 

(C.26) 



Just as equations (|C 13|) and (|C.14[) . equations (|C.21|) and (|C.23|) are integrated by two quadratures, defining 
by the way the constant of motion Q\: 



-40i 



(C.27) 



For Pi to be classically defined we need C 2 + Q\ > and C„ > or Q\ > 0. The domain of variation of e 2 ^ 1 is 



0, JC 2 + Qi+C 5 



e" 2 ^ G 



C 2 + Qi,C () + JC 2 + Q 1 



if Qi > , 

if Qi < and C g > . 



(C.28) 
(C.29) 



The first case corresponds to a bounce defining asymptotic Kasner's exponents; the second one to an oscillatory 
behaviour of the metric component. 

So we obtain the expression of the three variables 0i, 02 and 03 . Then the remaining three equations can 
be integrated by quadratures. For completeness we display them hereafter: 



The solutions: 

After elementary integrations, we obtain when K. 2 — S23 2 > : 

' K \p\ cosh[V/C? - s 23 ^ - t_)] + p S23 



03-02 = In 



B(t) := I J e^-M dt 



sgn(p) arctan 



K 2 - S 23 2 



/ K, - sgn(p) g 23 
K + sgn(p) s 23 



(C.30) 



ta,nh[x/lC 2 - s 23 2 (t - t-)/2] 



and when JC 2 — s 2 3 2 < and p S23 < 



ft = In 



/C |p| cos[Vs 23 2 -/C^(i - +PS23 



/C 2 - S23 2 



(C.31) 



23(i) = sgn(p) arctan 



' s 2 3 - sgn(p) /C 



tan[Vs 23 2 -/C 2 (i-i_)/2] 



s 23 + sgn(p) /C 

The function Bit) furnishes the time evolution of the non constant spin tensor components 

S [i2] = L cos[B(i)] , 
5 [i3] = Lsm[B{t)}. 



(C.32) 
(C.33) 
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Let us notice that, by fixing appropriately the initial values of the spinorial mode, we also have : p e^ 2 /33 ' +S23 = 
/Csin[2B(i)]. 

In the same way we obtain the expression of /3%. When Qi > it reads as: 

I [" ^JCj + Qi cosh[VQi"(i -*.)]- C g 



when Q\ < and thus C g > \Qi\, we obtain : 



Combining these equations with eqs (|C.26j) and eqs (jC.30l or IC.31[) we obtain the explicit expression of (3 2 and 
(3 3 . The functions v\ 2 and v\ 3 are then obtained by integration of eqs (|C.24p , IC.25|) where the spin tensor 
components are expressed thanks eqs (IC.321 IC33|) . The various integration constants p , . . . , Q, L are not all 
independent. The Hamiltonian constraint imposes: 

k 2 - pC 2 - - s 23 2 - L 2 = k 2 - K 2 - L 2 , (C.36) 

which is the classical analogs of the mass shell condition (|5.82[) , Obviously there are different qualitative 
behaviors of the metric components, according to the relative values and signs of some integration constants. 
From Eq. (|C.26[) and eqs (|C.30[[C".34[) we immediately obtain Kasner exponents (see next section). The function 
B(t), obtained from Eq. (IC.30[) . varies of a finite amount, and the spin tensor components S^ 12 ! and M 13 ' rotate 
by a finite angle. But contrary to what happens for gravity coupled to bosonic fields, the non-diagonal metric 
coefficients v\ 2 and v\ 3 , which are driven by these spin tensor components blow up exponentially. 

The solutions given by eqs (|C.31|) and (|C.35|) are of a different nature. Here the spin tensor components S 1 ' 12 ' 
and i!?' 13 ! oscillate. The metric component exp[— 2/3i] also oscillate between two positive values. The coefficients 
exp[— 2^2] and exp[— 2^3] are given by the product of a non-vanishing oscillatory factor and the exponential 
exp[2fci]. Accordingly they blow up when kt — > +00 and collapse to a singularity when kt — > —00. The 
evolution of the non-diagonal terms are different according to the subspace they concern. The coefficient ^23 
follows the derivative of P3 — [3 2 , and remains bounded, oscillating. The coefficients vyi and ^13 grow up like 
an exponential, but modulated by an oscillating factor. Notice that this last solution absolutely requires a non 
vanishing spinor field. 



(C.34) 



c 2 a 



Qi cos[V=3T(t-t.)] 



(C.35) 



Collision rules 



The collision rules, i.e. the transformation law of the Kasner exponents of the metric, when they exist (i.e. 
for Qi > 0), are directly read from the asymptotic behaviour of the solutions displayed in the previous sec- 
tion. These solutions allow one to discuss both collisions on a gravitational wall or on a symmetry wall. The 
arbitrariness of the constants f* and i_ reflects the independence of the order of the collision processes. 
The kinetic matrix is, with respect to the velocities {/3i, (3 2 , ^3, i>\ 2 , ^13, ^23 ]v 
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(C.37) 



The Hamiltonian constraint implies that asymptotically the trajectories in the /3-space are non-spacelikc: 



(3-f3 = \{w 2 



TaJ 2 + , — -{Wl + TU 2 

8 



™ 3 ) 2 = - 1 -(L 2 + S23 2 ). 



(C.38) 



4G 



with a mass term in agreement with [23 i. In terms of the velocities {A}, the collision rules are, for a collision on 
a gravitational wall {0i, ^2,^3} ^ {—0i, 02+2 01, 03+2 0i} and on asymmetry wall 01,02,03} l— ^ /?3, ^2}- 
In terms of momenta we obtain : 

{wi, -072, nj 3 } i-> {-mi + 2(w 2 + m 3 ), w 2 , m 3 } and {mi, ZU2, 1^3} i-> {^i, 073,^2} (C.39) 

These mappings conserve the value of the " mass" (|C.38[) . As is well known, the timclike nature of the trajectories 
makes that, after a finite number of collisions, the (3 "billiard ball" will end up on a worldline which does not 
catch up anymore the receding walls. 



The Dirac sector 



It remains to solve the Dirac equation (|A.29[) . taking into account the expression of the metric and spin com- 
ponents obtained. It is given by: 

X = N(t){ X , U} = \ 7 12 + 7 iS + (S m + e fe ~* n7 23 ) ^ - e" 2 ^) X (G.40) 

This classical equation for the rotation of the spinor \ is more complicated to solve than the quantum problem 
posed by the quantum type II Hamiltonian (|5.76p because it represents the rotation of the body-frame object 
X w.r.t. the space frame. This rotation is the combination of several different precessions and rotations 
corresponding to the various terms in the equation above. We can simplify a bit the study of this combined 
rotation by "encoding" the last term (linked to the interaction with the gravitational wall) in a change of spinor 
variable. Indeed, using the commutation of 75 with the spin matrices 7^ , its effect is taken into account by 
replacing, for Q\ > 0, x by the variable with: 



X = exp 



75 



ip = exp 



1 



arctan 



V 



ci 



c„ 



•tanh[v / S7(i-t*)/2] 75 



4>, (C.41) 



X = exp 



;| lJC 2 g + Ql+C g 



— arctan 
2 



tan(V^eT (t - U)/2) 75 



(C.42) 



when Qi < 0. 

Note that, after a collision on the gravitational wall, the spinor undergoes a total rotation around the 
"direction" 75 = 707 123 (associated to the gravitational wall) given by the matrix : 



8g := exp 



7r/2 - i arctan(yQi/C 9 ) ) 75 



(C.43) 



In particular, in the limit where vQi 3> C g , we get a rotation by an angle of j around 75 = 707 , in agreement 
with the sharp (gravitational) wall limit studied |22) . 

Let us finally discuss the effects of the other terms in the rotational evolution of the spinor, which, in terms 
of the object ip, are contained in the equation: 

jp = i (V i2 l 7 i2 + S [i3] 7 13 + (S [23] + e fe - fe m 23 ) 7 23 ) (C.44) 
For arbitrary values of L, it seems difficult to express simply the spinorial mode but assuming L — we obtain: 

m = (cos[i S23 * + \ B(t)] + sin[i S23 * + \ B(t)] 7 23 ) </>o • (C.45) 

This result exhibits a rotation in spinor space which is the sum of a uniform, continuous precession (with angular 
velocity S2 3 /4), and of the rotation B(t) (which spans a finite angle in the infinite time t € (—00, +00)). If 
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S23 = things further simplify. [Note that the two conditions L = and S23 = cancel the rotation of the body 
frame w.r.t. the space frame, and leave only the effects of the two dynamical walls: the gravitational wall 2/3i 
(studied above), and the symmetry wall /?3 — j3% that we consider next.] We obtain Bit) = arctan{exp[/C(i — t-)]} 
and after the collision the spinor will be rotated as : 

ip[t = -00] := O s tp[t = -00] . (C.46) 

in accordance with the collision law on a symmetry wall obtained in Ref. [22j . 

Thus to conclude this section let us remark that the bosonic part of the m = and = Einstein-Dirac 
system, under the assumption of a cosmological model of Bianchi type II, is completely integrable, in terms 
of elementary functions. For the fermionic part we have obtained explicit solutions only when simplifying the 
problem by, essentially, restricting the continuous rotation between the space frame and the body frame (to 
which the spinor is attached). However, as mentioned earlier it is not necessary to know precisely the spinor 
dynamics to elucidate the gravitational one. As discussed in the main text, the quantum dynamics of this 
system is simpler and can be solved in an exact manner. 



exp 



-^1 
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